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Abstract. Can high energy physics be simulated by low-energy, non-relativistic, 
many-body systems, such as ultracold atoms? Such ultracold atomic systems lack 
the type of symmetries and dynamical properties of high energy physics models: 
in particular, they manifest neither local gauge invariance nor Lorentz invariance, 
which are crucial properties of the quantum field theories which are the building 
blocks of the standard model of elementary particles. 

However, it turns out, surprisingly, that there are ways to configure atomic 
system to manifest both local gauge invariance and Lorentz invariance. In 
particular, local gauge invariance can arise either as an effective, low energy, 
symmetry, or as an "exact" symmetry, following from the conservation laws in 
atomic interactions. Hence, one could hope that such quantum simulators may 
lead to new type of (table-top) experiments, that shall be used to study various 
QCD phenomena, as the confinement of dynamical quarks, phase transitions, and 
other effects, which are inaccessible using the currently known computational 
methods. 

In this report, we review the Hamiltonian formulation of lattice gauge theories, 
and then describe our recent progress in constructing quantum simulation of 
Abelian and non-Abelian lattice gauge theories in 1 + 1 and 2 + 1 dimensions 
using ultracold atoms in optical lattices. 
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1. Introduction: Quantum Simulation of High Energy Physics 

1.1. Quantum Simulations 

Quantum simulation [T] is a relatively new physical field, concerning the simulation 
of quantum systems with other quantum systems. The original idea goes back to 
Richard Feynman [2j, who argued that since nature is quantum, it should be simulated 
by quantum computers. However, it was "left aside" for about two decades, as 
no candidates for playing the role "quantum simulators" were available. Following 
Feynman’s pioneering work [2], it has been shown that universal quantum computers 
Hi could in principle simulate general many-body systems [Sj as well as quantum 
fields [Bj, serving as digital quantum simulators. 

Thanks to the enormous progress carried out over the last decades on the road 
towards "Universal Quantum Computers", several quantum systems are now highly 
and precisely controllable in the laboratory, mainly belonging to atomic physics 
and quantum optics. These are, for example, ultracold atoms in optical lattices 
0 [HI EE m HI] i including Bose-Einstein condensates; Trapped ions [HUISUIllinHIB]; 
photonic systems EH, Rydberg Atoms m and more. Thus, although a full-fledged, 
universal quantum computer has not been realized yet, one may utilize these systems in 
order to create analog quantum simulators, which are indeed non-universal simulators, 
but nevertheless, allow for quantum simulation of quantum systems. This is done by 
mapping the degrees of freedom of the simulated systems into those of the controllable 
simulating system. Such quantum simulators provide simple realizations for quantum 
models and phenomena which may be, otherwise, inaccessible, either theoretically, 
experimentally, or numerically, amendable for direct observation of many-body physics 
phenomena. Using such simulators, one could study, experimentally, the spectrum and 
dynamics of the simulated system, and therefore these simulators may be regarded as 
analog quantum computers, specifically tailored to the simulated systems. Quantum 
simulation also has the possibility of realizing physical models which are "unreal" - 
not believed to be found in nature. 

One could also propose and realize digital quantum simulators, which fulfill the 
same tasks, but using quantum computation methods of qubits and quantum gates. 
This work, however, deals with analog simulations. 

Thus, quantum simulation is currently a rapidly growing physical area, involving 
multidisciplinary research, both theoretical and experimental. The theoretical research 
focuses on mapping the simulated physics into the quantum simulating systems, 
considering possible realizations and required approximations of the simulated model, 
as well as methods to control and get the most out of the simulating systems. Hence the 
multidisciplinarity: it involves both the simulated and the simulating physics, which 
may belong to two totally different physical disciplines, at least from a traditional 
perspective. The experimental research focuses on improving the capabilities and 
controllability of the simulating systems, also by developing new techniques and 
simulating systems, and of course, on realizing quantum simulations, serving as table- 
top experiments, with great success. 

Having been realized using AMO physics, models of Condensed Matter and 
statistical mechanics are the natural candidates for quantum simulations. These 
may be Hubbard models, spin chains and Heisenberg models and others - only few 

examples are [7] [TTJ] [2UJ [2TJ [2211231 [Ml [Ml (Ml HU [Ml [Ml (3HI • Effects of artificial, 

external gauge potentials, experienced by cold atoms, have also been considered 
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for quantum simulation, among more topological effects, such as the Aharonov- 
Bohm effect, Berry’s phase, topological insulators and the quantum Hall effect 
[511 155115511551 [5111551 [5511571 [5511551H51 pm 115] . A significant example to the power 

of quantum simulation in condensed matter physics is the experiment for simulating 
the phase transition in Hubbard models m ■ 

However, quantum simulations are not limited only to condensed matter physics, 
and they involve also other physical areas, such as gravity, including Black-Hole 
Hawking [15] and Unruh [11] radiation [IS] 05] [IT] [15] 05] [55] [5TJ [55] [55] . Relativistic 
physics has been considered as well, including, for example, Dirac Equation (including 
Zitterbewegung and the Klein paradox) [51; ;55, 55, 57], the famous Fermi two-atom 
problem [55; 55 > neutrino oscillations m and the Majorana equation [611 [55] . The 
next obvious step is to consider quantum field theories as well. 

Quantum simulations of quantum field theories involve these of vacuum 
entanglement of a scalar field [65, [54] as well as the interacting scalar-fermionic 
theories of Thirring and Gross-Neveu [65J. Fermions in lattice gauge theories have 
been considered as well, for example in the quantum simulations of Axions and 
Wilson fermions [55], Dirac fields in curved spacetime EH, general simulations of 
QFT with topological insulators [55, and Schwinger’s pair-creation mechanism [55] . 
Simulations of nontrivial geometries and topologies and extra dimensions, using 
internal degrees of freedom [75] up have been suggested too, as well as simulations of 
some supersymmetric theories, such as [75|[75||7I]. 

1.2. Quantum Simulations of High Energy Physics 

Another class of quantum simulations is those of high energy physics (HEP), which 
involve dynamical gauge fields. 

Gauge fields are in the core of the standard model of particle physics [75] [75]. 
These fields, through their local symmetry, induce the interactions among matter 
particles, which are fermions - or, phrased differently, the excitations of the gauge 
fields are the so called gauge bosons, force carriers, or interaction mediators. (In the 
current report, we shall not discuss the Brout-Englert-Higgs mechanism 132 El], which 
introduces a scalar field, not playing the role of a gauge field.) 

Each gauge theory is based on local gauge invariance - invariance to local gauge 
transformations, generated by the so called gauge group. Such groups may be either 
Abelian or non-Abelian: for example, QED (Quantum Electrodynamics) is an Abelian 
gauge theory, based on the group 17(1), while the strong interactions are described by 
the theory of QCD (Quantum Chromodynamics), which has the non-Abelian local 
gauge symmetry of SU( 3). These symmetries induce the interactions among matter 
particles, which have different gauge charges - the electric charge in QED , or the non- 
Abelian color charge in QCD, for example. In the case of non-Abelian groups, the 
gauge bosons (e.g., QCD's gluons) also carry charge of their own, unlike the Abelian 
photons of QED. 

The standard model of particle physics involves a local gauge symmetry of 
SU (3) x SU (2) x U (1), corresponding to the strong (QCD), weak and electromagnetic 
(QED) interactions respectively. All the gauge groups in this case are continuous 
groups - Yang-Mills theories m- 

The study of gauge theories, within and beyond the framework of the standard 
model, employs many theoretical techniques. These include, of course, the 
perturbative approach of Feynman diagrams [55] 151] [55] 155] [55], which has been 
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applied with great success to quantum electrodynamics. However, within QCD , such 
perturbative methods face a problem. In high energies, or short distances (such as the 
scale of deep inelastic scattering) perturbative techniques work fairly well, since, due 
to asymptotic freedom |85| (and the earlier results of the parton model and Bjorken 
scaling [MUSH), in this regime the quarks behave essentially as free, pointlike particles 
(within the hadrons). On the other hand, at low energies, or large distances, the quarks 
are subject to the so-called quark confinement [88j which forbids the existence of free 
quarks and rather binds them together into composite particles - hadrons. This limit 
is nonperturbative, and other techniques must be applied in order to study it. 

The most powerful approach that has been developed for non-perturbative QCD 
effects is lattice gauge theory eeiiibiiii. It allows probing and calculating many 
important quantities and phenomena of the theory, both analytically, as the lattice 
suggests a regularization method, and numerically, mostly using Monte-Carlo classical 
simulation methods. Such studies have led to significant and remarkable results, such 
as the low-energy spectrum of QCD and the CKM matrix elements m, several 
results concerning the quark-gluon plasma fS2J |jJ5] S3, [SJ[ and the deconfinement 
phase transition at finite temperature [96] 97.; ;98j - just to name a few. 

Monte Carlo calculations, nevertheless, are limited in certain important case. 
For example - the computationally hard sign problem [99, that limits calculations in 
regimes with a finite chemical potential for the fermions - which is problematic, for 
example, for the study of the QCD phases of color-superconductivity or quark-gluon 
plasma [931 1100 , 101J , or the fact that the Monte-Carlo calculations are of Euclidean 
correlation functions, and not of real-time dynamics. Indeed a full study of quark 
confinement with dynamical charges, and other nonperturbative phenomena in QCD 
in 3+1 dimensions is still lacking. 

Another successful avenue for the exploration of such theories has been the 
development of "toy models" - simpler physical theories which capture the essential 
physics in question. Such studies include, for example, confinement in Abelian theories 
in various dimensions (including on the lattice) |1021110311104111051 15511551 1106111071 
HTJ51 HOE or QCD 2 , the 1+1 dimensional version of Quantum Chromodynamics, 
which has been a target of interest over the last decades, as a natural playground 
for nonperturbative calculations, and as a basis for analytical understanding of real- 
world, 3+1 QCD ( QCD 4 ). There, several nonperturbative methods have been proven 
fruitful and useful for derivation of the spectrum of SU(N C ) theories. First, in the 
large N c limit, firstly discussed (and solved) by’t Hooft ffl~01 1111] . where the mesonic 
spectrum of the theory is revealed, manifesting confinement of quarks. Other methods 
to study the hadronic spectrum of QCD 2 have been DLCQ (discretized light-cone 
quantization) [ 112| and current algebras [ 113] . Witten’s non-Abelian bosonization 
m has helped to gain insight on the strong coupling Baryonic spectrum and quark 
content |ll5t 11161 1117 1 i 118 1. The relation of the fermion mass with confinement in 
these theories has been studied, for example, in [1191 1120] . 

Thus, new methods of calculations in gauge theories, in particular QCD , 
could be of a great help. Quantum simulation methods may be candidates for 
that. Once a quantum simulator for a gauge theory is built, one could use 
it to observe, experimentally, the otherwise inaccessible physics of these special, 
nonperturbative regimes. Such quantum simulators may be the quantum computers of 
HEP calculations: allowing for real-time dynamics and including real fermions rather 
than Grassman variables, suffering from the sign-problem, they suggest a way to 
overcome the problems of numerical lattice gauge theories mentioned above. Besides 
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that, reconstructing the standard model (in this case, out of atomic building blocks), 
or any of its sub-models, could help in understanding the basic interactions and 
symmetries of the theory. 

Quantum simulation of high energy physics, involving gauge field(s), may be both 
theoretically and experimentally harder than simulating other physical phenomena, in 
condensed matter physics, for example. In fact, three basic requirements K2U must 
be met in order to obtain a quantum simulator of a gauge theory: the simulation 
must include both fermionic and bosonic degrees of freedom (matter and gauge 
fields/particles); It must have a Relativistic, Lorentz invariance (a causal structure); 
and it must involve local gauge invariance , in order to obtain conservation of charges, 
and of course - the required interactions. These requirements are not trivially met, 
however, they must be satisfied, as without them, the theories in question could not be 
achieved. As we show in this paper, these requirements may be satisfied and fulfilled 
in certain configurations of atomic systems. 


1.3. Proposals for Quantum Simulations of High Energy Physics 


So far, only theoretical proposals have been suggested for quantum simulations of high 
energy physics. 

The first group of works contains condensed matter systems in which gauge 
invariance emerges, which may serve as quantum simulators for Abelian gauge theories. 
These include, for example, an effective emergence of a 17(1) spin liquid in Pyrochlore 
1 122) : A ring exchange model using molecular states in optical lattices, yielding a 
Coulomb phase 17(1) in the limit of no hopping [ 123] ; The effective emergence of 
artificial photons in dipolar bosons in an optical lattice m; a digital quantum 
simulation of spin-| 17(1) gauge theory, as a low energy theory of Rydberg atoms 
m-, and the emergence of gauge fields out of time reversal symmetry breaking for 

5 


spin-§ in honeycomb optical lattices 

Quantum simulation, using ultracold atoms in an optical lattice, of continuous 
QED has been suggested in | T2'6] . However, most of the works so far concentrate 
on quantum simulations of lattice gauge theories, which are more suited for quantum 
simulation using optical lattices. 


Simulations of pure-gauge compact QED (cQED , see section 2.2.1 ) include a 


simulation of the Kogut-Susskind Hamiltonian in 2+1 dimensions using Bose-Einstein 
condensates in optical lattices m, as well as a simulation of a truncated, spin-gauge 
model, using single atoms in optical lattices |128| . 

Simulations of U( 1) models with matter have been suggested as well - either 
for an Abelian link model 12!) 130 . 131| - a 1+1 dimensional simulation of the 
lattice link Schwinger model [132] , or a generalization of the spin-gauge approach, to 
include dynamical fermions in 2+1 dimensions cm The above proposals for cQED 
simulations with ultracold atoms introduce gauge invariance as an effective symmetry, 
by constraining Gauss’s law in the Hamiltonian, with a large energy penalty. Another 
approach for simulation, utilizing fundamental symmetries of the ultracold atoms, is 
proposed in m, for 1+1 and 2+1 dimensional cQED , with or without dynamical 
matter. 

Recent proposals for such simulations also utilize other systems - such as 
superconducting quantum circuits [13411135| and trapped ions |136| , both simulating 
the Abelian link model. A digital simulation of a Z7(l) gauge theory is presented in 
cm using Rydberg atoms. 
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Simulations of discrete lattice gauge theories have been suggested as well, 
including discrete Abelian (Zjv) and non-Abelian (D/v) models using arrays of 
Josephson junctions [I38j . a digital simulation of Z 2 with Rydberg atoms |137| and 
an analog simulation of Zjv, with an explicit construction for Z3 in an optical lattice 
m- A way to truncate the Hilbert space of a U( 1) theory in a manner similar to 
Ztv was proposed in m- 

Quantum simulations of Non-Abelian gauge theories with continuous groups 
(SU( 2) Yang Mills models, with possible generalizations to U(N), SU(N)) have 
been proposed too - Analog simulations with cold atoms utilizing prepotentials 
[1301 114111142 1 1143 1 or rishons in the link model [1441 1145 ] have been suggested in 
|14fi| and respectively, where the first utilizes hyperfine angular momentum 

conservation to obtain local gauge invariance, and the latter starts with a larger 
symmetry group, which must be broken. A digital simulation of an SU{2) gauge 
magnet |l30 . 148j has been suggested in [ 149] , while a digital quantum simulation of 
an SU (2) link model with triangular plaquettes, using superconducting qubits, has 
been suggested in |M]. A new way for imposing gauge invariance using dissipation 
(Zeno effect), with an explicit construction for a non-Abelian model, has been given 

in nsu. 

In general, the dimensions of the analog quantum simulations with optical lattices 
are restricted only by experimental and technological considerations. Thus, the above 
suggestions for 2 + 1 dimensional systems, using optical lattices, should apply, in 
principle, also to d + 1 systems, given, of course, that the required technological 
challenge is met. The 1 + 1 quantum simulations mentioned above which use optical 
lattices may be generalized to further dimensions but, mostly, without the magnetic 
(plaquette) terms, and thus are simply quantum simulations the strong limit of 
lattice gauge theories (see section [ 2 ] for an explanation of these terms). For other 
simulating systems and/or digital simulations, the case may be different, due to 
possible dimension-dependent properties of the simulating system. 

For a comprehensive review on the recent progress in quantum simulation 
proposals of link-Rishon models, the reader should refer to jl52[ . 

Throughout this work, h = c = 1, and the Einstein summation convention is 
assumed. 

2. Lattice Gauge Theory: A brief review 

Lattice gauge theories (5911901 [T53 , T541 3~5Bl 1156 are formulations of gauge theories on 
a discretized space, or spacetime. They were originally invented by Wilson, as a tool for 
the study of the quark confinement problem [55] . It is an important tool in high energy 
physics, and especially in QCD (Quantum Chromodynamics - the theory of the strong 
interactions) as it allows performing nonperturbative calculations (numerically, using 
Monte-Carlo methods) and thus provides insights into the perturbatively-inaccessible 
regions of QCD, due to the running coupling. 

Besides addressing the confinement of quarks (or charges in analogous Abelian 
theories, as will be later discussed), lattice gauge theories have enabled the 
computation of other important quantities in QCD, such as the hadronic structure 
and spectrum, as may be read in the reviews [9011101119T] and many others. 

In the following review of lattice gauge theory, we follow the Hamiltonian 
(canonical), Kogut-Susskind formulation of lattice gauge theories [ 15711591190 1 (rather 
than the more "conventional" Euclidean approach). This is since the Hamiltonian 
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Figure 1 . The lattice structure. The gauge degrees of freedom occupy the links, 
whereas the fermionic matter - the vertices. 


language is much more natural in the context of atomic and optical physics, used for 
the simulating systems. 

2.1. Hamiltonian Formulation 

2.1.1. Hopping fermions and global gauge symmetry. Consider fermions on a d 
dimensional spatial lattice, with nearest-neighbor hopping. The Hamiltonian takes 
the form 

H = (ipkp n+ k + h.C.) (1) 

n n ,k 

where n £ 7L d label the lattice’s vertices and k £ {1,..., d} labels the lattice’s directions; 
k is a unit vector (where the lattice spacing is a = 1) in the fcth direction. The 
Hamiltonian consists of local “mass” terms and hopping terms (see figure |T|) . 

?/> n are spinors of some unitary gauge group G, in some representation r [|] For 
example, f/’ n are merely fermionic annihilation operators in the case G = (7(1). For 
G = SU(N ), in the fundamental representation, these are iV-component spinors, 
containing N annihilation operators: 


Ipn 


( V’n.l \ 

Ipn,2 


V V’n.JV / 


( 2 ) 


where in Hamiltonian Q, and in what follows, summation on the group indices is 
assumed, e.g. 

V’nV’n = (3) 

i 

The Hamiltonian (|T]) is invariant under global gauge transformations of the gauge 
group G. That means, that if we pick some group element V G G (which is unitary) 


| Note that we only consider here the gauge degrees of freedom; Each component of these gauge 
spinors may also be a spinor in terms of the spin or flavor degrees of freedom, but we disregard it 
here. 
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and perform the transformation 

V'n » Vlpn ; Ipi *• (4) 

the Hamiltonian is left intact [§] This corresponds to a conservation of the total number 
of fermions, 

Ntot = (5) 

n 

which is obvious, as all the terms of the Hamiltonian annihilate one fermion and create 
another one. 

2.1.2. Local Gauge Symmetry. Next, let us lift the symmetry to be local, by 
generalizing the transformation to 

-> VnV'n ! V’n > ( 6 ) 

- i.e., each vertex is assigned locally a unitary element of the gauge group. The mass 
terms of ([l]) are local and hence remain intact, but the hopping terms vary under the 
transformation, for example 


^nV’n+k — 5 ► i’lVlV n+ tV» n+k (7) 

The gauge symmetry will be restored if we introduce a connection t/ nj fc, located 
on the link emanating from the vertex n in the fcth direction. is an element of 

the gauge group in the appropriate representation (corresponding to the spinors to 
which it is coupled), a unitary operator [jj] undergoing the gauge transformation as 

u a , k —> V a U n , k Vl +t ( 8 ) 

then, if the Hamiltonian is changed to 

HgM = (^i^n.fc^n+k + h ' C ) ( 9 ) 

n n,k 

the gauge invariance is, indeed, restored, and the gauge transformation leaves it intact. 
The name Hqm stands for the gauge-matter interactions described by it. 

We denote the local gauge transformation generators by G l n . Due to the gauge 
symmetry, they commute with the Hamiltonian, 

[H,G l n \= 0 (10) 

and thus the physical Hilbert space TL is also gauge-invariant, i.e. divided into sectors 
of eigenvalues of G l n (see figure [ 2 ]) : 

•h=©-h(R}) (11) 

{Qn } 

§ The group element V in Q is, of course, written in the appropriate representation (the same as 
the spinors). However, note that the equation is representation independent. When a representation 
is chosen it must be identical for the spinor and the group element. 

|| The unitarity has many advantages and is used in the context of conventional, Wilsonian minimally- 
coupled lattice gauge theories. However, the transformation (symmetry) properties may be satisfied 
if the operator is not unitary as well, as we do in our simulations 112811133111211 [1461 . 
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• ■ ■ 

! n= © n{{< 

£}) 


Figure 2. The physical Hilbert space of a gauge invariant theory is divided 
into separate sectors, corresponding to different static charge configurations GD- 
represented by boxes. The gauge invariant dynamics does not mix among them 
\13\ . These are two consequences of gauge invariance ©■ 


Each such sector TL ({<z^}) has a set of eigenvalues {(?„}, called static charges , such 
that for every \ip ({q l n })) € TL ({<&}), 

Gj 1 |^(R}))=4|V’({4})) (12) 

Note that as the generators do not commute for a non-Abelian group (this shall be 
clear in the following subsection, where they are explicitly defined). The static charges, 
thus, can only be defined by the eigenvalues of a subset of some commuting operators. 
The above equation is the generalized Gauss’s Law, and the dynamics does not mix 
these sectors, i.e. 

| H\ij) ({g^})) « ( 13 ) 


2.1.3. Dynamical Gauge Fields. We shall identify the connections U £ with the 


gauge group elements, i.e. introduce the group parameters 9 l n k - Hermitian operators, 

such that in a given representation r, with the matrix representation of the group 

( r ) 

generators T> , 


U (r) k = 

n, k 


AGl 


p( r ) 


(14) 


where a summation on the group indices i is implicit. We shall assume next 
working in the fundamental representation (generalizations to other representations 
are straightforward), and disregard the representation indices in this case: U nt k = 

One could then define operators which are conjugate to the group parameters, 
to serve as generalized momenta. In case of the Abelian group U (1), where = 
e *0n,f% one could define the conjugate momentum to be proportional simply to 0 n ^. 
But generally, when G may be non-Abelian, one must pay attention of the non¬ 
commutativity of the group in doing that. Thus, we could reformulate the Abelian 
case and define the conjugate momentum there as 


Ln,k — iGn.kG^ k — k Un,k 


(15) 


The differential representation would be 


L 


n,fc — 


. d 

1 d6 nk 


(16) 
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from which the canonical commutation relation 


[0n,fc) L n ,k] = i (17) 

results, and from which, as it is the generator of angular translations, 

[L n ,k,U n ,k\ = U nk (18) 


Thus a reasonable choice for the dynamic part of the gauge field would be 

Hi CX J2 L lk 
n ,k 


(19) 


in agreement with the Abelian version of the Kogut-Susskind Hamiltonian (see below). 

Let us now generalize it for non-Abelian groups |158j . Due to the non-Abelian 
nature of the group, one could define two sets of such operators, left and right - 

L n ,k = —iU n ,kU^ k ; Rn,k = —'i-UikU-n.k (20) 

(for an Abelian group L = R, but this does not hold in general). One can expand 
these operators in terms of the group’s representation matrices, 


L n _k — L l „ uTi 


k^i i 


Rn.k = Rl k-Ti 


( 21 ) 


and then construct the differential forms of the operators , |i?“ fc |. These 

operators are called left and right generators of the group, since they satisfy the 
group’s algebra U] 

( 22 ) 


[Li k ,U ntk ] = T x U n , k 

[K,k,u a ,k] =u n , k T, 

T'n ,kiL J n k = ~~ 

[K , k ,Rl,k\=ifniL l n , k 


r i 

-^n.fei-n-n ,k 


= 0 


(23) 

(24) 

(25) 

(26) 


where {7)} are the group’s generators and f\ji are the group’s structure constants p~| 
Generators and group elements from different links commute, of course. 


Note that due to the definition (201, one could prove that 


Rn,k — U i.L^kUn.k 


(27) 


and thus if we calculate the trace in group space we find out that 


tv(Rl k )=tv(Ll tk ) (28) 

Note that the right and left operators may be defined in other conventions as well, resulting in 
different signs in the group’s algebra. 

f In a general representation r, these equations would be generalized in a straightforward manner, 

[ (r) 1 (r) (r) (r) 

L l n k , k \ = J U ^ k , where Tr stands for the matrix representation of the generator in r , 

(r) 

and k stands for the matrix representation of the group element. If one wishes to explicitly write 
the group’s indices, the equations will look like, for example, fc , (i7 n) fc)^J = (TiU n ,k) 
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This can be reformulated using the generators, 

E«fc) 2 = E( L -,fc) 2 = L n,fc (29) 

i i 

giving us the electric part of the non-Abelian Kogut-Susskind Hamiltonian, 

H e cx ]Tl4 lfc (30) 

n,k 

We can refer to jl/W|, j | as the left and right electric fields, conjugate (in some 

sense) to the vector potential. Their difference along a link corresponds to the charge 
carried by it: no charge in the Abelian case, and the color charge in the SU ( N ) case. 
The gauge transformation generators are then given by ina 

G; = E ( £ L - K-h.) ~ «n (31) 

k 


where Q “ is the dynamical charge, which will be introduced later. 
The Abelian version of the generator is merely 

Cn = ^2 — ^n-k,fc) ~ Q™ 

k 

= div n L n fc Qn 


(32) 


in which the discrete divergence is defined, and from which the relation to Gauss’s law 
is obvious. 

One would also like to consider a gauge-field self interaction. Gauge invariance 
forces us to consider only group traces of products of group elements along close 
paths. The smallest such interactions are the so-called plaquette terms, involving unit 
squares of the lattice, from which we obtain the magnetic part of the Kogut-Susskind 
Hamiltonian [157] 

H b oc 5Z (t r (UiUzUluS) +/i.c.) (33) 

plaquettes 

with the group elements usually in the fundamental representation (even if there are 
charges in other representation which requires the use of other representations when 
coupling them to the gauge field) which simplifies in the Abelian case to [59] 

Hb oc E cos (0i + d 2 ~ 0 3 - 0 4 ) (34) 

plaquettes 

where the numbers in both the equations are according to the plaquette convention 
presented in figure [3] Note that this agrees with the terms of the Wilson action [55], 
and thus with the “usual” procedure of obtaining the Hamiltonian (by a Legendre 
transformation of the Lagrangian). 

We have obtained the terms of the Kogut-Susskind Hamiltonian only up to some 
proportion factors. One can find them by taking the (classical) continuum limit. 
Doing that (or using a Legendre transformation of the Lagrangian) one obtains the 
Kogut-Susskind Hamiltonian 

H KS = \g 2 Y>l^Yg 2 E (tr (lhU 2 Utut)+h.c) (35) 

n,k plaquettes 

where g is the coupling constant. 
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Figure 3. The convention of numbering the links in plaquette interactions (see 
equations (|33|l , (|34|) ). 


2.2. Abelian Theories 

2.2.1. Compact QED (cQED) The first theory we shall consider is Compact 
Quantum Electrodynamics - cQED |89] . This theory’s continuum limit is Quantum 
Electrodynamics (QED). Its lattice version is compact, involving angular variables. 

This is an Abelian, continuous gauge theory, with the gauge group U (1). Thus, 
the group elements are simply unitary phases 

U n , k = (36) 

as we have already discussed, 8 n k are analogous to the vector-potential, and the 
conjugate momenta are T n ,fc = — i - - U( 1) (planar) angular momentum operators. 
The canonical commutation relation [0 n ,fc,-bn.fc] = i ensures that t/ nj fc is a raising 
operator of the electric field [L n ^,U nt k] = U n} k, and U^ k is the lowering operator, 
both are unitary. 

The Kogut-Susskind Hamiltonian takes the form 

= E cos(9 1 + e 2 -0 3 -e 4 ) (37) 

n ,k ^ plaquettes 

where the indices inside the plaquette term are with respect to the convention of 
figure [3] Also note that the argument of the cosine corresponds to the lattice curl, 
manifesting that the continuum limit of this term gives rise to the magnetic energy. 

Let us consider the local Hilbert space on each link. One could work in the basis 
of electric (flux) eigenstates, which is a good basis in the strong coupling limit g 2 1, 
as shall be discussed when we consider confinement. The eigenstates of L on a single 
link satisfy 

L | to) = m | to) (38) 

where to takes any integer value, either positive or negative. In coordinate (angle) 
representation, the wavefunctions are 

(9) = (0\m) = -E e ime (39) 

The Abelian Gauss’s law takes the form 

On = E (E.fc “ -^n-k.fc) “ Qn 

k 


div n L n ,fc Qn 


(40) 
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and thus, in this basis, the gauge-invariant states satisfy (without dynamical matter) 

Gn |0) = “ m n~k ,k) 10 ) = 10) (41) 

k 

where q n is the (C-number) static charge. 

What are the interactions? 

(i) The plaquette interaction, in this language, changes the electric flux on a 
plaquette, such that in two links the flux increases by a single unit, and in the 
other two it is lowered. This is done only in two possible orientations, which are 
the only possibilities which leave the eigenvalue of G n intact in all the four vertices 
of each plaquette - thus, this is the “physical interpretation” of the commutativity 
of the plaquette interactions with Gauss’s law. 

(ii) If we wish to introduce dynamical charges, staggered |1571 159| . for example]^] 
with a Hamiltonian of the form 

Hm = eY + H c ■) + m Y ^ 0n0n (42) 

n ,k n 

where {n^} are the indices of the vertex n. we may interpret the gauge invariant 
interactions as ones which raise the dynamical charge on one edge of the link, lower 
it on the other side, and raise/lower the electric flux on the link in a way that 
Gauss’s law, now including the dynamical charges, will hold in both the edges. In 
this staggered representation, even vertices represent “particles”, with a positive 
mass to, and odd ones - “holes”, with a negative mass —to. In fact, the negative 
vertices may be interpreted as the Dirac sea: suppose we shift the energy by a 
constant, and measure the mass of the odd vertices with respect to —to. Then 
an occupied odd vertex corresponds to a state with zero mass - no particle at 
all, and a vacant odd vertex corresponds to an anti-particle with mass to. This 
corresponds to the canonical transformation of second quantized Dirac fields, 
where the holes’ creation operators are replaced with anti-particles’ annihilation 
operators. We can also define the charges of the particles: the particle would 
have charge 1 (in units of the fundamental charge, of coupling constant g), and 
the anti-particle - charge -1. This means, in the language of fermionic operators, 
that 

, 1 / Vn fe \ 

Qn = ^ y 1 - M) k J (43) 

and this completes the particle-anti-particle picture. 

These two possible interactions may be understood as the meaning of local 
gauge invariance. Note that the only two possibilities of contracting gauge invariant 
operators are traces of products of group elements along closed loops, or along open 
loops, bounded by two fermionic operators on the edges. These two interactions 
are the shortest such interactions possible, and thus the most local ones. The same 
considerations apply as well to other lattice gauge theories. 

* Note that staggered fermions are not the only possible recipe for discretization of fermionic fields, 
required by the doubling problem 115911160119011153 ! ' 1541 1T55I I156 : . The different approaches, which 
are out of the scope of this paper, vary in both quantitative and qualitative details, and have their 
own advantages and disadvantages. Yet, when a quantum simulator is designed, one should consider 
which type of lattice fermions to use. We mostly use the staggered fermions formulation, but other 
approaches should be possible as well. 
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2.2.2. Zpf. Another Abelian gauge theory of relevance is Zjv • This is a discrete 
group, and thus this is not a Yang-Mills theory, and its structure is a little different. 
Nevertheless one can use this group to formulate a lattice gauge theory and observe 
interesting physics. This group is highly relevant for the confinement in QCD , as Z 3 is 
the center of this group, and the center is responsible for confinement (large-distance 
phenomena) [ 161] . 

On every link of the lattice, we define two unitary operators, P and Q [ 162| . such 


that 

P^P = QlQ = 1 

(44) 

which satisfy the Zjy 

algebra, 



P N = Q N = 1 ; ptgp = e iSQ 

(45) 


where 6 = 

Let us define the eigenstates of P as 

P\m) =e im5 \m) (46) 


there are N such states, 


to £ 


N —1 


N- 11 

~ 2 ~ J 


(47) 


for an odd N (the generalization for an even N is straightforward). Q is a unitary 
ladder operator, 

Q |to) = [to — 1) (48) 


with the cyclic property 


Q 


N- 1\ 
~ 2 ~ / 


N- 1\ 
~ 2 ~ / 


(49) 


alternatively, one may use the eigenstates of Q, and then P will be a unitary raising 
operator, cyclic as well. 

One can define Hermitian operators E, A such that 


P = 


e iSE 


Q = e z 


(50) 


and in the N —> 00 limit they will correspond to cQED 's conjugate electric field and 
vector potential. 

The Hamiltonian of this theory takes the form [11621 

H = (- P + P lk) -\Y, (Q 1 Q 2 Q 3 Q 4 + H.c .) (51) 

n,k plaq. 

with local terms and plaquette interactions, which tend to the Abelian Kogut-Susskind 
Hamiltonian ( |37[ ) as N —> 00 . The plaquette indexing convention is, again, according 
to figure [3] 

One can define static modular charges the vertices 

qn = e- iSm (52) 


and then the Gauss’s law is given by 


( 53 ) 
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where 

Gn = H P IU P > ( 54 ) 

1 + l- 

with l + being links which start from n (positive links) and /_ - ending there (negative 
ones). 

A description of gauge theories with finite groups (both Abelian and non-Abelian), 
with dynamical, fermionic matter as well, is found in [163J. 


2.3. Non-Abelian Yang-Mills Theories: SU(2) as an example 


As an example of a continuous Non-Abelian gauge theory, we examine next the SU(2) 
gauge theory ( 1571 ;90J. SU(2) (the rotation group) has three generators, with the 
structure constants 


fijk 

— ^ijk 

(55) 

right algebras 



— RijlL n k 

(56) 


(57) 


The pure-gauge Hamiltonian is (35). One may choose the representation of the 
U matrices, and let us do that in the fundamental representation of SU (2) (j = |). 
In this representation we work with the Pauli matrices, 



(58) 


and thus 


I'n.k — 6 2 


SI 


(59) 


- which are 2x2 matrices of operators. The reader may find an explicit expression for 
this operators, using the formulation [ 163| . when we discuss the quantum simulation 
of this model. 

The links’ local Hilbert space, in the flux basis, is now described using three 
quantum numbers |jmTO , ). These are eigenstates of the operators 


L 2 | jmm') = R 2 \jmm') = j (j + 1) \ jmm') 


(60) 


L z | jmm') = m | jmm') (61) 

R z | jmm') = m' \ jmm') (62) 

A mechanical interpretation of this Hilbert space may be a rigid body, described 
in both body and space coordinates. In both of them the total angular momentum is 
equal, but in these two frames of reference the components of the angular momentum 
will be measured differently, as they form two separate, commuting algebras - the body 
and the space algebras, corresponding to the two algebras we have here [ 1571 Il64 j. 

By acting with the group elements, which are rotation matrices, on the | JMM') 
states, an amount of angular momentum is given to the state (or taken from it). 
Unlike in the U( 1) case, where the angular momentum is described by a single quantum 
number, here the representation is also a dynamical quantity: for example, acting with 
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U n .k in the fundamental representation on a | JMM') state means adding a ^ quanta 
of angular momentum to the state, resulting in a superposition of J + \ and J ~ \ 
states (ordinary addition of angular momentum). Quantitatively this is expressed by 

CSS] 


U mm , | JMM') = C+ 


J + -, M + m, M' + m' ) + C_ 


J — M + m, M' + m' 


(63) 

where C± depend on the appropriate Clebsch-Gordan coefficients [1551 1166 1 . Thus the 
representation of the state, the quantum number J, is a dynamical observable here, 
unlike in continuum field theories, where each field has a constant, well-defined, static 
representation: this is due to the fact that lattice gauge fields are not Fock-space 
bosons [ 163] , 

If one wishes to include dynamical fermions as well, the interaction with which 
shall be described by a Hamiltonian of the form where ip n are group spinors. 
Thus, if we pick U to be in the fundamental representation, this will also be the case 
for the spinors, which will contain now two components - each Lorentz component will 
be such a spinor. With staggered fermions we shall have such a single color spinor on 
each vertex. In this case the charges will simply be (1571190] 


Q 


(Oij V’n.i 


(64) 


satisfying the group’s algebra. Note that the “Dirac sea” state can be obtained here as 
well (as a no-interactions vacuum). Again, we measure the mass of each particle with 
respect to — m, and a fully occupied vertex, due to the fermionic nature of the matter 
particles, corresponds to zero charge. Thus the Dirac sea is the case with empty even 
vertices and fully occupied (by two fermions) odd vertices [ 15?] . 

The description of other SU ( N) groups is similar; Changing the group means 
changing the local Hilbert spaces, but the Hamiltonian has the same form. A 
description of Hamiltonian lattice compact QCD {SU (3)) is given in [lOOj . 


2-4- Physical content of the models 

Gauge theories in general, including ones on the lattice, exhibit an interesting phase 
structure. These include, for example, QCD 's phases of color-superconductivity 
or quark-gluon plasma [Ml llOO i 1101] . whose treatment using the usual Monte- 
Carlo methods utilized in lattice gauge theory is rather problematic, due to the 
computationally hard sign problem m that limits calculations in regimes with a 
finite chemical potential for the fermions. 

Another interesting and important property of gauge theories is the existence of 
confining phases, dependent on the coupling constant and other parameters. We shall 
review quark confinement in pure gauge theories in detail below. 


2.4-1- Quark Confinement. Free quarks are not found in nature, but rather bind 
together and form hadrons, such as baryons and mesons. The reason to that is the 
physical phenomenon of quark confinement, manifested in QCD but also in other 
gauge theories: in large distances, or low energies, the static potential between quarks 
is linear in the distance, and thus it is energetically unfavorable (and practically 
impractical) to separate bound quarks. It is a highly intriguing concept, forming 
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a fundamental aspect in the understanding of the Hadronic spectrum, involving 
nontrivial low energy physics, which cannot be addressed perturbatively. 

Schwinger showed, by exactly solving 1+1 dimensional QED with massless 
charges mm, that vector mesons can be massive if massless fermions are screened 
by the vacuum polarization, obtaining massive photons. Further studies nnnnna 
have shown that only massive photons, no electrons, exist in the asymptotic states 
of this model. The latter work has also shown, in accordance with the expected 
confining behavior, that electrons arise in deep-inelastic-scattering scenarios - i.e., 
short distances and high energies, and they behave like free pointlike particles (agreeing 
with the “parton model”, or the asymptotic freedom of QCD, where the confined 
quarks are asymptotically free within the hadrons - at short ranges and high energies). 
The existence of massive photons is in accordance with the fact that in terms of 
statistical mechanics, confinement is a gapped, disordered phase, while the ordered 
phases do not confine |1681 1158 j . 

A significant step in the study of confinement was done by Wilson, who introduced 
lattice gauge theory for this purpose [55]. Using loops in lattice gauge theory, which 
were later called “Wilson Loops”, he showed that 3+1 dimensional compact QED has 
a confining phase (for strong coupling, g 2 1), which binds charges together and 

forbids the existence of free electrons and positrons, as well as a Coulomb phase (for 
weak coupling, g 1 <C 1), as in continuous QED. It was shown later (see subsection 


2.4.31 that 2+1 dimensional cQED confines for all values of the coupling constant 
m confinement in compact QED was identified as a topological effect, having to 
do with instantons and magnetic monopoles. This result was supported by other works 


as well HrmniMnnui . 

Similarly to the “charge confinement” in Abelian theories, confinement also takes 
place in non-Abelian gauge theories, such as non-Abelian Yang-Mills theories (as 
QCD) as well as the discrete theory of presented earlier. It was argued |ll6l i that 
large-distance phenomena, involving the confinement of quarks in QCD, is related 
with the group’s center - Z 3 , which makes the phase structure of Zjv interesting 
P21[TM1IT70] : The theory confines in the strong coupling limit (A — > oo). In 2+1 
and 3+1 dimensions it is not the only phase. In 3 + 1 dimensions, due to the self¬ 
duality of the Zjv theory, there is a phase transition from electric confinement (strong 
coupling) to magnetic confinement (weak coupling), for N < N c (N c ~ 6). For 
N > N c there is a third phase, with no confinement at all. In 2+1 dimensions there 
is a phase transition to a non-confining phase in the weak limit |1681 l('i2 . However, 
for N —> oo the theory shows the phase transition at g = 0 m, in accordance with 
the single confining phase structure of 1/(1) (the N —> oo limit of Zjv)- 

Both in 1/(1) and SU(N), one can show, using perturbation theory, the 
confinement of charges/quarks in the strong coupling limit, as will be done in 
subsection [2~j~2| The confinement of static charges is manifested by “electric flux 
tubes” connecting them, which get thicker as the coupling constant decreases. But, 
outside the strong limit, confinement (in cases it holds) becomes nontrivial and requires 
nonperturbative techniques, and for Yang-Mills theories it is still an open question: 
in fact, the proof and calculation of the Mass Gap in Yang-Mills theories, which is 
related to quark confinement, is one of Clay’s institute Millennium Prize problems. 
For example, in the cases confinement was proven, it was mostly done for static charges 
- and not dynamical fermions. 

Dynamical fermions shall enhance electric flux-tube breaking. If one stretches a 
flux tube, at a certain point it would be energetically favorable to create two pair of 
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“new” quarks, or charges, breaking the a long flux tube into two smaller ones. 


2.^.2. Strong Coupling Lattice Gauge Theories. An interesting limit of lattice Yang- 
Mills theories is the strong coupling limit g 2 1, in which the Kogut-Susskind 
Hamiltonian (351 may be treated perturbatively: the magnetic, plaquette part Hb 


may be treated as a perturbation to the electric He = Hq. [|] 

Consider two opposite static charges, placed on the vertices n and n + Rk of 
a lattice. What is the ground state of this static charge configuration in the strong 
coupling limit? In this limit, products of local (links) electric field states form a good 
basis, as they are eigenstates of He, and thus we obviously must find the electric field 
configuration with minimal energy, which satisfies the static Gauss’s law for these two 
charges. 

The no-charge vacuum of this limit satisfies 




> l)) = <g> |0) 

n,k 


(65) 


- on each link the state satisfies 


L 2 |0) = 0 


( 66 ) 


- no electric field at all. 

If we excite a link (using U) the zeroth order energy is increased by a unit of L 2 : 


L 2 U |0) = [L 2 ,t/] |0) 



(67) 


where T[ is the matrix representation of the group’s generator in representation r, 
and C <2 ( r ) = X) is just a C-number - the eigenvalue of the Casimir operator in 

i 

this representation. 

Thus, for example, in the case of cQED , we get [89 

L 2 U |0) = U |0) (68) 

and in the case of SU (2), for charges in the fundamental representation j = \ [95] . 


L 2 U |0) = 



(69) 


or in any other representation r, 

L 2 C/ r |0) = j r (j r + 1) U r |0) (70) 

This means that if we have N links with a single excitation, the energy of the 
state (in zeroth order) will be N 'j C 2 (r). In order to respect Gauss’s law (and fulfill 
gauge invariance), the eigenstates of He in this charge configuration will contain a 

jj It is worth to comment that strong coupling calculations were initially performed in the Euclidean 
formalism. The first to do that was Wilson m ■ 
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flux line of U |0) states, ranging between the charges. The minimal energy is obtained 
for the direct line, and thus the ground state is: 


n+(i?—l)x 

I Ro,mm>) = I I] I l fl > 


(71) 


El This flux-line is the so-called electric flux tube , manifesting confinement. 
The zeroth order energy of this state is 


Eo | No, mm' ) = He \Ro,mm ') = ypC 2 (f) R \Ro,mm ') (72) 

and thus we get that the strong coupling limit confines, with a string tension (the 
ratio between the static potential and the distance) of 

2 

<? strong = y C 2 (r) (73) 

Such a flux-tube connecting the appropriate pair of them is nothing but a 
meson. One could also create a long flux tube connecting static charges and watch it 
dynamically breaking into shorter, energy favorable tubes, by dynamic pair creation. 
Also, as mentioned before, outside the strong limit, but within the confining phase, 
the state shall be modified such that the flux-tube is thicker (wider than just a thin, 
single lattice line). 

2-4-3. The Weak Limit of cQED . As g decreases, perturbative corrections must 
be taken into account, until perturbation theory breaks down and nonperturbative 
methods of the weak regime must be utilized. As an example, we shall comment 
about the weak limit of cQED. 

In 3+1 dimensions, one can use continuum limit Wilson loop arguments (89j, as 
the weak limit corresponds to the continuum limit, and obtain the Coulomb potential. 

In 1+1 dimensions, on the other hand, the theory confines for all values of 
the coupling constant, as it is solvable and has no dependence on the coupling 
constant g. This was shown (for the massless case) in the continuum by Schwinger 
PHI HQ3), but similar results apply also in the lattice compact case in m (as 
well as in other, more recent, lattice studies, such as using tensor network states 

fTTSlITTSl H7il H75llT7BllT771ir75| l. 

In 2+1 dimensions interesting physics takes place, involving instantons and 
topological effects. Due to magnetic monopoles, there is a mass gap in the weak 
limit, corresponding to quark confinement. This was shown by Polyakov [ 106] using 
Euclidean partition functions, by Drell, Quinn, Svetitsky and Weinstein, eng using 
variational Hamiltonian approach, and by Banks, Myerson and Kogut jl07j . using the 
Villain approximation. The three approaches were compared by Ben-Menahem [109 
who also demonstrated that compactness is necessary for confinement in lattice QED. 

ft The indices in the state are another manifestation of the fact we disregard the charge degrees of 
freedom in this context. If one introduces the matter operators n , the operator creating the state out 

/ n+(H— l)x \ n-\-(R— l)x 

of the full vacuum would be J2 V+m II U n , x V'n+Rx.m' = W. II t+s+n+Rx 

mm' V n=n / W 

which we easily identify as a gauge invariant string. 
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In some sense, the confinement mechanism of 2+1 dimensional cQED is more 
related to the confinement of quarks in QCD, as both of them are topological effects. 

The discussion above applies for zero temperature. For a finite, nonzero 
temperature, a Coulomb phase arises for weak values of the coupling constant, with a 
phase transition to the Coulomb phase at some value of T, g mm- 

3. Ultracold Atoms in Optical Lattices : A brief review 

3.1. Generation of an Optical Potential 

Optical lattices are created by “external”, classical lasers, and virtual processes to an 
excited atomic state dunum]. 

Consider a two level atom (disregard the other levels), whose levels are |a) - some 
meta-stable level (in which the atom is initially prepared), and |e) - an excited one 
(see figure |4|. We may write its Hamiltonian in the form 

p 2 

Hatom = x-1+ |e) (e| (74) 

Zm 

The atoms are usually alkaline, therefore they have a single valence electron, which 
simplifies their treatment by many aspects. 

The atom experiences an external, classical laser field, of the form 

E (x, t) = E (x) e~ iLJt e (75) 

Switching to a system rotating with the laser’s frequency w, and defining the detuning 
S = w e — w (see figure [4]) , we get 

Hatom = £-+6 |e)(e| (76) 

Zm 

The interaction between the laser and the atom is given in the dipole 
approximation, which is valid if we assume that the electric field’s amplitude varies 
slowly compared to the atomic size (in space) and 1/w (in time; Note that here, 
for simplicity, the amplitude was not even considered as having a time dependence, 
but one may generalize to a time-dependent amplitude, of the form E (x, t)). This 
interaction has the form 

Hdipoie = -d E (x, t) (77) 

where d is the dipole moment of the operator; After transforming to the rotating 
frame, performing a rotating wave approximation (neglecting the fast-rotating terms) 
and plugging dij = e ■ (i| d \j) , the dipole moment matrix elements (projected on the 
laser’s polarization), one obtains the interaction Hamiltonian 

Hdipoie = -d ea E (x) \e) (a| + h.c. = | e ) ( a | + h.c. (78) 

where 

n (x) = -2 E (x) d ea (79) 

is the Rabi frequency ( d ea = d ae , as the dipole operator is Hermitian). 

We assume that the detuning S is sufficiently large - the coupling is non-resonant, 
and that transitions between the atomic levels are practically impossible. However, 
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Figure 4. The energy levels of the trapped atom. |o;) is the meta-stable level, for 
which the optical potential is effectively obtained by eliminating the excited level 
|e). The levels are separated by frequency cj e , and they are coupled non-resonantly 
with a laser of frequency u the detuning is 5 = u e — uj. 


virtual second-order transitions to the excited level and back are possible, and thus 
we get the effective Hamiltonian (for the |a) subspace) 



(80) 


where 


K p (x) 


l»(x)| 2 

46 


(81) 


is called the optical trapping potential, and it is an effective potential experienced 
by the atom in the |a) state. In fact, this is merely the AC Stark effect. Note that 
we have actually “traced out” the internal levels, due to the energy restriction (large 
detuning) and what we now have is a configuration space Hamiltonian. Choosing 
the effective field properly, one could generate an optical lattice: a periodic trapping 
potential, allowing the trapped atoms to fill its minima. Since V op (x) ex |i£(x)| 2 , a 
standing wave of the form E (x) = E 0 cos (kx) would result, for example, in an optical 
potential of the form V op (x) = Vo cos 2 (kx); One could add lasers in more than one 
direction to obtain a two- and three-dimensional optical lattices of various geometries. 

Next we shall consider an optical lattice with many atoms. This shall be done in 
the context of second quantization, for which we are interested in the single-particle 
eigenfunctions of the trapping Hamiltonian. 

Due to the lattice symmetry, the single particle wavefunctions are merely 
Bloch waves. These states are not localized, therefore if one wishes to work in 
the number basis of local minima (Fock basis), a transformation to local single 
particle wavefunctions is required. This is fulfilled by the orthonormal basis of 
Wannier functions, w n (x — xo), which are the discrete Fourier transform of the Bloch 
wavefunctions Uq' 1 (x) over the lattice mmm 


w n (x - x 0 ) cx lq ' x °uW (x) 

q 


(82) 


These states are labeled by n, which is the energy band, and are each spatially centered 
around a different local minimum Xo of the optical potential. These are the desired 
local site (minimum) wavefunctions. 

Working in low enough temperatures, it is reasonable to include only the lowest 
energy band, wq (x — xq). 
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One can expand the second-quantization “wavefunction” of the condensates in 
terms of the localized Wannier functions <f> n (x) = wq (x — x n ) and single particle 
annihilation operators a n (bosonic or fermionic, depending on the statistics of the 
atoms), both defined at the site (minimum of the potential) n, as 


$ (x) = ^a n </> n (x) 


(83) 


and construct from which a usual second quantization Hamiltonian density: 

U = (x) + V op (x) + V T (x)^ $ (x)+ J d 3 x' J d 3 x& (x') & (x) V (x - x') 4> (x) $ (x') 

( 81 ) 

where Vp (x) is an external trapping potential, varying slowly compared to the optical 
lattice potential V op (x), and V (x — x') is some two-body interaction, which will be 
specified later. 


3.2. Multi-species Lattices and Spinor Condensates 

It is also possible to utilize more atomic levels and transitions and create a “species- 
dependent lattice”: i.e., with different minima for different atomic species, which may 
be, for example, different atomic hyperfine levels [T7inrT80l 1261181 ITT]. 

The atomic hyperfine structure is given by an SU(2 ) Hilbert space, characterized 
by two quantum numbers F,m F . Just like “regular” angular momentum, 

F 2 \F,m F )=F{F + l)\F,m F ) (85) 

and 

F z |F, m F ) = m F |F, m F ) (86) 

F is defined as 

F = I + L + S (87) 

where I is the nuclear spin, and L, S stand for the valence electron’s orbital angular 
momentum and spin respectively (recall that we are dealing with alkaline atoms, with 
a single valence electron), and thus the Lie algebra 


[Li ■ Fj] iCijkFjz 


( 88 ) 


is satisfied. Each species |a) can then experience a different optical potential (x), 
and these potentials form a superlattice. 

Such lattices / condensates are sometimes referred to as spinor lattices / 
condensates. 

As an example, consider s 'Rb and 23 Na have a nuclear spin I = |. Thus, in 
the S-wave hyperfine manifold (which is relevant for ultracold atoms, as explained in 
subsection 3.31, F = 1 or 2. 

We denote the internal levels by Greek indices, and the second-quantization 
expansion generalizes to 

(89) 


(x) — ^ (^x) 


for each species, where the Wannier functions may be species dependent or 
independent - depending on the lattice structure and symmetry and the choice of 
these levels. 
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The Hamiltonian density is now 

72 


n= (x) ( 5 «n f _ Yl + v « (x) + (x)') + ( x )) ( X ) 

a,/3 ^ ' 

j (Px'&a (x') (x) V a p lS (X - X') $ 7 (x) (x') 


(90) 


ot, / 3 , 7 , <5 


where fl Q/3 (x) corresponds to a Rabi coupling of atomic levels by more lasers. 


3.3. Scattering of Ultracold Atoms 


How cold is ultracold? The ultracold regime is defined as the one in which the lowest 
partial wave governs the scattering. It is a dilute and cold regime. 

We are dealing with the so called weak interaction regime , involving the ground 
state of a gas of either bosonic or fermionic atoms. For ideal bosons, we assume 
that they all occupy the lowest energy band (Bose Einstein Condensate). A phase 
transition from Bose-Einstein condensate takes place in a finite temperature (in three 


dimensions), when the thermal de-Broglie wavelength A t = \JjrNp (in H = 1 units) 

approaches n -1 / 3 - the average inter-particle distance (n is the number density). For 
fermions, the states are full up to the Fermi energy, ep: for temperatures satisfying 
ksT < ep the gas is considered degenerate |9|. 

When the energy is low enough - i.e., the de-Broglie wavelength of the atoms is 
much larger than the effective range of the potential in, the interaction potential in 
(841 is accurately described by the pseudopotential: 


V (x - x') 


2770 ^( 3 ) 

m 


(x-x') = f S& (x-x') 


(91) 


For bosons the scattering is governed by the S-wave channel. For identical fermions, 
on the other hand, due to Pauli’s exclusion principle, there is no S-wave scattering JS]; 
This is not a problem for the scattering of two different fermions (non-identical). 

The density of the atoms, n, is typically 10 12 — 10 15 cm -3 . This gives a typical 
interparticle distance n -1 / 3 of the order of 0.1 — 10/rm, and the scattering length is 
in the range of few nanometers [S|. The scattering length of 81 Rb, for example, has 
on order of magnitude of 100 times the Bohr radius, varying slightly according to 
the different internal states ms]. The gas is sufficiently dilute for ?r 1,/3 a <C 1, in 
case of bosonic repulsive interactions (a > 0), as Vna 3 is the small parameter of the 
Bogolyubov theory, describing weakly interacting Bose gases [5]. 

We thus obtain a simplified version for the single-species case Hamiltonian density: 


n = & (x) + V op (x) + P T (x) j $ (x) + |$t (x) (x) $ ( x ) $ (x) (92) 

after plugging in the Wannier functions and creation/annihilation operators (831 one 
may integrate over the Wannier functions. Using the overlap integrals 


e n = J d 3 X(j)* n (x) + V op (x) + V T (x)^ (f n (x) 

Jmn = J d 3 X<t>* m (x) + V op (x) + V T (x)^ (f n (x) 


(93) 

(94) 






CONTENTS 


26 


Umnki = 9 J d 3 x(j)* m (x) (/>* (x) <^ fc (x) <j)i (x) (95) 

where m,n, k,l label the lattice sites. The most general Hamiltonian, taking into 
account all possible overlap integrals, is 

di = Y. J ^nala n T 'y ) (96) 

m,n m,n,k,l 

3-4- Multiple Species Scattering 

In case of multiple species, the most general form of the Hamiltonian is 

H= E E C C&>L,a4,/J«fc,7“^(97) 

ra,n,a:,/3 m,n,k,l,a,(3,‘y,5 

however, if we remember that the species are different hyperfine states, and take 
into account the hyperfine rotational symmetry (conservation of the hyperfine angular 
momentum in atomic collisions) we can simplify this Hamiltonian further, as we shall 

do next H7S1IISD1. 

Suppose we wish to consider the scattering of two atoms with distinct F,m F 
values: for example, let us scatter an Fi atom on an F 2 atom. As a result, we shall 
get a pair of two atoms with their own values. However, taking into account 

the rules of addition of angular momenta, one should be aware that the scattering 
may take place “through” an “intermediate state” with several values of total angular 
momentum Ft- Thus, the effective scattering potential will take the form 

9_- 

y a ,p,y,s (x - x') = =-d (3) (x - x') '^2a FT (PF T ) atPni s (98) 

Ft 

where we sum over the possible values of total angular momentum Ft, and consider, 
for each value, its corresponding scattering length a Fr ■ Pf t is the projection operator 
onto the subspace of total angular momentum F T . In case of different atomic masses, 
the appropriate reduced masses should replace the mass m. 

The projection operators can be represented using functions of Fi -F 2 . Physically, 
that may be understood as a consequence of rotational invariance: terms which depend 
on the F values of two atoms and conserve the hyperfine angular momentum must be 
constructed out of powers of these terms. Mathematically, 

Fi • F 2 = \ (f! - Fl - F1) (99) 

and as F 3 , F| are practically constants for our purposes, a function of Fi ■ F 2 is in 
fact a function of F^, and thus one can construct the projection operators out of this 
scalar product. If there are n possible values of total hyperfine angular momentum, 
the projection operators are polynomials of the form 

n— 1 

P Ft = {Fl ■ F 2 ) fc (100) 

k =0 

and eventually obtain 

W(x-x') = « (3) (x-x')f((F 1 -F 2 )‘) 

Z \ / Q!,p,7,0 


( 101 ) 
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where the {gk} coefficients depend simply on the scattering lengths {a Fj .} and 
the coefficients in the projection operators’ expansion, {Gp Tt k}- involved. The 
scattering lengths are tunable using Feshbach resonances, either magnetic or optical 

jTHTl ITH^l [THTH fTHTTl [T71 [H?T1 

We can also express the projection operators, and thus the scattering potential, 
in terms of the Clebsch-Gordan coefficients: 

V a4 3, 7 ,5 (x - x') = S (3) (x - x') 5> r (Fi, m F1 = a; F 2 , m F2 = /3| F T , M F ) x 

Ft 

(Ft, M f | Fi,m F i = 7 ; F 2 , m F2 = 8) 

which explicitly manifests, as a consequence of hyperfine angular momentum 
conservation, that a + /3 = 7 + 5 , and the {C Ft } coefficients depend on the scattering 
lengths. 

As an example [179 . ISO] , one may consider a spinor condensate of F = 1. The 
possible values of total angular momentum in collisions are naively Ft = 0 , 1 , 2 , 
but since the atoms are bosons, and the states must be symmetric, we remain with 
Ft = 0,2. Thus the scattering lengths we need are ao and a 2 , and the highest power 
of Fi • F 2 we need is 1. The projection operators are 

P 0 = ^(l-F 1 -F 2 ) (102) 

and 

P 2 = \ (Fi • F 2 + 2 ) 

One can define 

47T 

9o = 7 ^ (2 a 2 + ao) 

and 

47 r 

52 = 3^ (a2_ao) 

and obtain the scattering part of the Hamiltonian density 

|E $ «( X )^( X )^( X ) $ “( X ) + f ( X ) F «7 $ 7 W) • (H ( X ) F/3«$5 (x)) 

a,/3 a, ,0,7, <5 

(106) 

where F are the F = 1 representation matrices of SU(2). 

3.5. Summary of the Available Interactions and Possibilities 

To conclude, in optical lattices one could use either bosonic or fermionic ultracold 
atoms, of species {a}, with single particle local wavefunctions at sites n , 4> n ,a (x) and 
ipn.a (x) for bosons and fermions respectively, and annihilation operators a njQ and c n<a 
for bosons and fermions respectively, satisfying the canonical commutation relations 



a n,a, a m ,P 

— ^nm^a/3 

(107) 

1 

Cn,on C m,/3 

* — $nm$a.(3 

(108) 


— {Cn,ot. •> ^m,^} — [^n,Q!5 — 0 

(109) 


(103) 

(104) 

(105) 
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Out of these operators and wave functions, the second-quantization held operators 
are constructed, 


(x) — ^ ^Qn,o:0n,Q; (^) 
n 

(110) 

^ oc (x) — ^ ^n,Qi^n,Q (^) 

(in) 


n 


and they are the ingredients of the atomic Hamiltonian, consisting of three main parts: 

(i) The single-particle terms, 

H 0 = [ d 3 * ( x ) Hoj'ka (x) + (x) H 0jb $ a (x)) (112) 

where Hoj,Ho,b are the fermionic and bosonic single-particle Hamiltonians, 
including the kinetic energy and the optical and trapping potentials. They depend 
on (or, from another perspective, affect) the shape of the lattice (through the 
optical potential) and thus the single-particle wavefunctions, 4> n ,a ( x ) and V’n.a ( x ) 

(ii) The scattering terms, 

Hsc = \ 9a0 7 S / ( X ) ^ ( X ) ^7 ( X ) ^<5 ( x ) ( 113 ) 

q:,/3,7,<5 

+ \ [ d 3 x^i (x) (j)\ (x) <£ 7 (x) (j> s (x) 

a,/3,7,<5 

+ \ 9a(hS J d 3 x^l (x) 4>l (x) <F 7 (x) (f)g (x) 

a,/3,7,<5 

which govern the two-body S-wave scattering (atomic collisions), and depend 
on the scattering coefficients g. They depend on the scattering lengths, and 
of course on the nature of collisions - fermion-fermion (FF), boson-boson (BB) 
or boson-fermion (BF), and the atomic species involved. These coefficients are 
not independent of each other, thanks to the conservation of hyperfine angular 
momentum, and are tunable using Feshbach resonances. 

(iii) Rabi (laser) terms, 


Hr = I d 3 x^ a (x) ^ (x) 

a,(3 

+ 

a,/3 

which are governed by external lasers and may induce desired hopping processes. 

One could also include transitions to and from molecular states, but these are not of 
relevance for this report. 

Current experimental technologies allow creating complex and controllable optical 
lattices. The range of parameters is wide, as exemplified by the typical numbers 
presented above: the ratio Vq/Er is tunable. Optical lattices may be created in 
various geometries and dimensionalities, using holographic techniques, for example. 


d 3 x<t >',L (x) fa (x) 
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The lattices may be very large (for example, 150000 lattice sites in [20]) and contain 
many atoms. The atoms may be distributed across the lattice in various ways: lattice 
sites may occupy few atoms, even single ones iaa Esa H2B ebb esj, or much more 
(for example, 80 in |192| . using tube-shaped sites). 

Common measurements are the so-called “time-of-flight” measurements, in which 
the lattice is switched off (suddenly or adiabatically) and the atoms are released, in a 
way that their spatial density is proportional to their momentum distribution, which 
is measurable using standard imaging methods [5]. Contemporary experiments also 
demonstrate with great success single-addressability of lattice sites and atoms m 
eehu ma eeh], which opens the way for various possibilities, both in measurements 
and initial state preparation. 

Further detailed information on the topic may be found in EH- 

4. A framework for Quantum Simulation of lattice gauge theories with 
ultracold atoms in optical lattices 

We shall review several methods for quantum simulation - i.e., physical realization 
using ultracold atoms in optical lattices, of three different lattice gauge theories: 
compact QED, SU(N) and Zjv- 

In order to simulate a gauge theory, three requirements must be fulfilled E2H: 

(i) The theory must contain both fermions and bosons. The fermions will represent 
matter fields and the bosons - gauge fields. This makes ultracold atoms natural 
candidates, as bosonic and fermionic atoms are practically free resources in these 
systems. Ultracold atoms in optical lattices are currently highly controllable 
systems. Of course, the nature of lattice fermions should be chosen carefully, based 
on both the experimental considerations and the desired theoretical approach for 
lattice fermions, as discussed above. 

(ii) The theory must be Lorentz invariant - i.e., to have a causal structure, as a 
field theory. This may seem problematic, but one can simulate the lattice gauge 
theories described above, that obtain this symmetry in the continuum limit. This 
goes along quite well with the structure of the simulating systems, which are 
optical lattices. 

(iii) The theory must have local gauge invariance , which is the symmetry responsible 
for gauge-matter interactions. This seems the most problematic demand, and the 
biggest challenge, as this does not seem to be a fundamental symmetry in such 
atomic systems, and has been one of the most important research steps. 

Thus, after considering these three requirements and recognizing that ultracold 
atoms in optical lattices should enable their fulfillment, such systems have been chosen 
as the simulating systems - building blocks for optical realizations of lattice gauge 
theories, opening the way to quantum simulations of lattice gauge theories. The 
main idea is to obtain the gauge invariant interactions which are not fundamental for 
ultracold atoms. We present two methods for that - one in which gauge invariance 
is imposed as a constraint - the effective method, in which the gauge invariance is 
only an emerging, low-energy sector symmetry, and the other one, in which gauge 
invariance is mapped into a fundamental symmetry of the atoms - conservation of 
hyperfine angular momentum in atomic collisions. 

In all the simulation approaches presented below, regardless of the specific 
implementation and representation of the gauge degrees of freedom, they are 
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Figure 5. (a) The lattice structure. The gauge degrees of freedom occupy the 

links, whereas the fermionic matter - the vertices. This structure is kept in the 
atomic quantum simulator - whose schematic plot is presented in (b): the circles 
and ellipses represent the minima of the optical potentials - Bosonic (B, yellow) 
on the links, and Fermionic (F, blue) on the vertices. 

(a) [H.GyXO (b) [/f,G‘ n ]=0 





Figure 6. Effective gauge invariance. The original atomic Hamiltonian (a) mixes 
among the Hilbert space sectors of different static charge configurations, as it is not 
gauge invariant: the generators of gauge transformations ( |31[) d o not commute 
with the Hamiltonian. Using the Gauss’s law constraint ( |115| >, one effectively 
eliminates the transitions among these sectors, resulting in the desired gauge 
invariant interactions, and the restoration of gauge invariance {To} by constructing 
an effective Hamiltonian to such a sector of interest. 


represented by bosons populating the links of an optical lattice. Fermions, in proposals 
involving dynamical matter, are always located in minima coinciding with the vertices 
(see figure]^)). 

5. Effective gauge invariance emerging at low energies 

First we shall describe a class of models, in which the gauge symmetry is not 
fundamental, but rather arises in a low energy sector, as an emerging, effective 
symmetry. In this method, the atomic Hamiltonian is not gauge invariant, but rather 
contains a large penalty term, demanding, as a low energy constraint, that Gauss’s 
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law will be satisfied. The constraint Hamiltonian takes the form 


He = A 


( 115 ) 


n 


where A is the largest energy scale in the system. This automatically divides the 
physical Hilbert space into sectors, corresponding to different eigenvalues of G n 
(311 - different static charge configurations, but the Hamiltonian interactions involve 
transitions among them. These are practically impossible, due to energy conservation. 
Thus, an effective Hamiltonian ; 193 T94| for the ground sector has been introduced in 
any of these proposals, giving rise to a gauge-invariant effective theory, which includes 
the desired interaction (see figure [6]). 

5.1. cQED in 2 + 1 dimensions - link interactions 

It is shown in m, as a proof of principle, that a quantum simulation of the Kogut- 
Susskind Hamiltonian for cQED is, indeed, possible, and this was the first work 
arguing that. This first proposal discussed only the pure-gauge theory, suggesting a 
method to observe confinement of static charges, by probing electric flux tubes. 

In this work, every link of a two-dimensional spatial lattice (labeled by the vertex 
from each it emanates (m, n ) and its direction k ) is occupied by a Bose-Einstein 
condensate, with a mean number of atoms Nq (uniform all over the lattice). We 
expand the local number operators AU n , 



(116) 


where Nq is a C-number and 5^ n - an operator, measuring the excess of the atomic 
population over (or below) Nq , which can take negative values as well - and represents 
the electric field. These operators are conjugate to the local condensates’ phases, 6 ^ n , 
playing the role of the compact vector potential. 

In a limit where the expectation value, as well as the uncertainty of <5^ n are 
much smaller than Nq, a quantum-rotor approximation may be introduced, allowing 
the expansion of local annihilation (or creation) operators as 


VN~ 0 e l8 


(117) 


a 


Neighboring links are populated by condensates of different atomic species, and 
thus direct, gauge-variant tunneling is eliminated. This is possible by the use 
of holographic masks, for example, for the creation of the optical lattice. Local 
scattering processes are responsible to the electric Hamiltonian, as well as, along with 
neighboring links scattering, Gauss’s law constraint. Using external lasers, desired 
hopping processes between neighboring links (diagonally) are introduced (see figure 
[7|. These terms are gauge-variant as well, but they induce effective gauge invariance: 
these hopping processes correspond to violation of the constraint, but second-order 
perturbation theory ties pairs of them together, forming the plaquette interactions 
(see figure [8]) . 

Experimental considerations (homogeneity of the scattering lengths) result in a 
sum-Gauss’s law, rather than with a divergence, i.e., 



( 118 ) 
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Figure 7. The lattice structure of m- Bosons of four different species 
(1-4) occupy the links such that neighboring links are occupied by different 
species, and thus nearest-neighbor hopping is avoided. Instead, diagonal hopping 
processes are induced by external lasers (denoted by black arrows on the right): 
1 «-»■ 3,3 •<->■ 2,2 «-»■ 4,4 •<->• 1. Scattering processes take place among the four 
species meeting in (and overlapping on) a vertex, and must be tailored properly. 


and thus a (theoretical) canonical transformation, involving inversion of signs, must 
be carried out in order to obtain the “right” cQED degrees of freedom, 



m+n q 



E k m , n = (“I) 


Being theoretical, this transformation should be carried out on the results of the 
measurements: a strong coupling electric flux-tube will thus be manifest by a line of 
alternating values of S!^ n between the charges. 

This proposal also allows leaving the perturbative (in terms of the simulated 
theory) strong-coupling regime, but to a limited extent, since the atomic excess 
must be small comparing to the mean population, which cannot be too large due to 
experimental consequences. However, in 2+1 dimensions, where confinement holds 
also in the weak limit, the simulation of the weak limit should work - at least 
qualitatively, even if not describing the exact weak limit of the full Kogut-Susskind 
Hamiltonian. This can be understood from the weak limit analysis of m- 

5.2. Truncated cQED model - link interactions 

The second proposal [1281 was a first step in the direction of a simpler implementation. 
Instead of Bose-Einstein condensates, each link is occupied by a single atom, belonging 
to each of 2£ + 1 possible atomic species, representing 2£ + 1 possible values of the 
electric field on a link: Here we have introduced the spin-gauge model, in which the 
electric field is truncated (taking values between — i and £), but the gauge symmetry 
is unaffected, which should lead, at least qualitatively, to similar effects those of as 
the Kogut-Susskind Hamiltonian. The U (1) angular momenta are replaced by SU (2) 
angular momenta in an integer representation £ with 2£+l levels, with the L operators 
of the electric field replaced by the truncated L z , and the unitary raising/lowering 
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Figure 8. Generation of plaquette interactions as effective second-order terms, 
as in m- First, an atom hops from a link to a neighboring one (diagonally). 
Then the Gauss constraint is violated, and a second process, restoring 

the constraint, must take place. This may either be another hopping process, 
with other two links, closing a plaquette, or the inverse hopping which results in 
“canceling” the first one, resulting in a constant term in the approach of EH (and 
in other gauge-invariant terms in the spin-gauge approach m, where the link 
operators are not unitary). Each plaquette is constructed out of four different 
second order processes. Note that the arrows correspond to the “modified Gauss’s 
law” |ll8[ ), with opposite signs on half of the links - i.e. to the language of the 
simulating systems (atoms) and not to the “right signs” of cQED , obtained after 
the canonical transformation t |119[ |. 


phase operators replaced by the non-unitary SU (2) ladder operators L±: 


U 


y/e(e + i) 


L+ 


( 120 ) 


in a confining phase, the electric field is supposed to take small values to, and thus 
with a large enough £, as well as the m <C l of confinement, the unitarity of the ladder 
operators is approximately restored, at least in the sense that when one acts with 
them on the | £m) states, no prefactor is obtained: 


1 


zL± | Im) = 




vW+ff 

In this representation, one obtains the Spin-Gauge Hamiltonian |128| . 
1 


Hsg = -rrT.-Tn.fc ^ 


2 g 2 P (i+1 f 


£ ( 

plaquettes 


-£'+,n.l-£'+ i n+l,2-k — : n+2,T—v 


( 121 ) 


H.c. 


( 122 ) 

We have compared this model to the Kogut-Susskind Hamiltonian (37). First, 
we have shown the more obvious result, that in the strong limit, as the states are 
eigenstates of electric flux (see subsection 2.4.2), one may use a truncated theory 
which suffices to describe the state of the system in the presence of static charges. 
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Conditions on £, depending on the magnitude of the static charges, have been discussed 
using perturbation theory considerations. It was shown what is the highest order in 
strong-limit perturbation theory, in which the spin-gauge states will correspond to the 
Kogut-Susskind ones, for a given charge configuration. 

Furthermore, we have addressed the nouperturbative, nontrivial weak limit as 
well, and compared numerically the single-plaquette results for a truncated Kogut- 
Susskind theory as well as the spin gauge theory to the analytical results of [108 . The 
results converge quantitatively to the analytical ones very quickly - for small values of 
£ (already for 2,3). 

An implementation for £ = 1 - three bosonic internal levels - has been introduced. 
It involves a three-scale hierarchy - the first, strongest constraint is responsible to 
the effective "angular-momentum" interactions [ 1951 and the second one to gauge 
invariance, in a similar manner to mi- Working with l = 1 limits the possibility of 
going into the weak limit. Nevertheless, the rapid convergence of the results in the 
weak limit to the ones of the Kogut-Susskind model, do show that implementations 
of the spin-gauge model for higher values of i shall serve as an accurate weak-limit 
simulation. Confinement and flux tubes take place here in the form of alternating 
atomic species along the line connecting the static charges - in the strong limit. 

Due to the non-unitarity of the L± operators, the effective constructions involves 
some undesired terms. However, they are gauge-invariant as well (as they satisfy 
Gauss’s law constraint) and their effect is discussed in |128| . 

Finally, in H33l> we have introduced dynamical fermions as well. Still considering 
cQED , in the spin-gauge model presented in paper [128] . but regardless of the 
implementation the gauge degrees of freedom, we have introduced naive fermions into 
the system, in 2+1 dimensions. 

Working with naive fermions in two spatial dimensions, each vertex contains two 
fermionic species, simulating a 2-component Dirac spinor, 


V’n = 


(123) 


with the local charge 


Qn = V’nV’n - 1 = N% + ~ 1 


(124) 


(similarly to the charge in the second quantized Dirac field. There, for the lower entry 
of the spinor, we swap the fermionic creation and annihilation operators, and thus 
obtain that this charge transforms to the well known — N^ form). 


The simulated Hamiltonian (besides the pure-gauge part, (122)) in 2+1 
dimensions, is 


H„ 


/Jf +1) 52 - H - C •) + (125) 

V V ' ) n.fc n,k 


with the Dirac matrices being Pauli matrices, aq : 2 = o x , y and f) = a z . The fermionic 
mass terms manifest explicitly that the particles and antiparticles are not divided 
among vertices as in the staggered formulation. 

Gauge invariance is still effective here. Inclusion of dynamical charges into the 
Gauss’s law constraints made it harder to implement, requiring also alternating signs 
to the dynamical charges. A detailed prescription for the implementation, based on 
fermion-boson scattering, was given, enabling the generation of this required constraint 
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- the interaction between the gauge field and the matter arises as an effective second 
order term as well: in order to satisfy the constraint, fermions can only hop along 
with a change of the electric field on the link. 

Generalization of this method to other types of lattice fermions is straightforward. 
One simply has to change the fermionic ingredients and tailor the desired hopping. 
The effective construction will "make sure" that only the gauge-invariant interactions 
are obtained. 

5.3. Effective Gauge Invariance - some general remarks 

The three proposals | 1'2 7 . 12811133] involved gauge invariance as an effective symmetry. 
This allows introducing gauge invariance to the simulating systems, but it is not the 
only possible approach, as shown in the subsequent papers (121 1 146 1. 

The Gauss’s law constraint involves scattering of atoms, representing the gauge 
field degrees of freedom, on several neighboring links. All the interactions must be 
with the same amplitude, which is hard to implement. This may lead to practical 
impossibilities to satisfy the Gauss’s law, for example, or to the unintended inclusion 
of fractional charges. 

The robustness of this method was considered by Kasamatsu et. al. fTMlIWj . 
who showed, mainly based on [T27] but also on other works of us and the other groups, 
that one may relax the fine-tuning demands of the simulating Hamiltonian parameters, 
and still get a gauge-invariant theory, involving Higgs fields. 

Another disadvantage of this approach applies to simulation of non-Abelian gauge 
theories. There, the generators of gauge transformations, required for the Gauss’s 
law constraint, are much more complicated (several non-commuting generators ©), 
which makes the constraint Hamiltonian flTT5t very difficult to realize. 

6. Local gauge invariance arising from many body interaction symmetries 

The basic idea of the second simulation approach we describe hereby, is that one could 
use local interactions of the many body system - atomic collisions - to achieve local 
gauge invariance. In particular, we have utilized the conservation of the hyperfine 
angular momentum in atomic collisions. This enables to obtain, as fundamental 
ingredients of the Hamiltonian, the on-link gauge-matter (boson-fermion) interactions 
of without any constraints or perturbation theory |121 ; : 146| . As a second step, 
in 2 + 1 (and more) dimensional setups, one may introduce the loop method |121| 
presented below, which allows the construction of the plaquette interactions effectively, 
but out of the already gauge invariant on-link interactions. 

The basic, on-link interactions ([£]) are obtained from boson-fermion collisions, 
utilizing the overlap of four Wannier functions on a link (two bosonic on the links 
and two fermionic on its edges). The only allowed scattering channels are those which 
conserve the hyperfine angular momentum. Thus, one can specifically choose the 
hyperfine levels playing the roles of the different simulated bosons and fermions, such 
that the F conserving processes will be mapped to the gauge invariant ones. This 
gives us the desired interactions, along with other ones, gauge invariant as well, which 
may be tolerated or tamed. 
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6.1. Compact QED - on-link interactions 

We shall show this procedure in detail for compact QED in 1+1 dimensions - the 
lattice Schwinger model. Besides its advantage as a "toy model" which allows a clear 
and simple presentation of the method, the quantum simulation of 1+1 dimensional 
QED with matter, in the manner presented below, has several advantages: 

(i) As gauge invariance is fundamental here, no use of effective interactions is required 
- i.e., the simulating Hamiltonian produces the simulated theory without any use 
of perturbation theory. This implies that the energy scales are larger, and thus 
the simulation should be faster and more robust to noise and decoherence. 

(ii) It is a possibly simpler implementation (compared with | j 127! and [128] ). 

(iii) The 1+1 dimensional Schwinger model is exactly solvable H021HD5I . Lattice 
results are available as well f [171 and much more, including new calculations 
involving matrix product states cm Along with 3, this make this model a 
possibly preferred model for a first physical realization, as one has analytical and 
numerical results to compare with. 

While the fermions are simply fermionic atoms in this simulation (staggered ones, 
to be specific, but they may be generalized to other lattice fermion representations), 
the bosonic fields are constructed out of the so-called Schwinger bosons, utilizing 
Schwinger’s representation of SU(2 ) [ 1981 Il99j : In this cQED simulation, each link 
may be populated by two different bosonic species, a and b, from which the algebra is 
constructed in the following way: 


tfb) 

; t = - (a? a + b^b) 

(126) 

a)b ; 

L_ = frl a 

(127) 


using this representation, along with the conservation of hyperfine angular momentum, 
simulations of the spin-gauge Hamiltonian may be obtained, but much simply than 
in the effective approach. Assuming that the bosons do not interact among links, £ 
is fixed and governed by the number of on-link bosons. If we consider a single link, 
whose left edge may be occupied only by the fermionic species c , and its right edge 
by d, we wish to obtain the interactions 

c l L + d + d) L_c = c^a^bd + d^b^ac (128) 

and this is possible by selecting 

mF ( a ) + tof ( c ) = mj? (b) + mj? (d) (129) 

(see figure [9]). Even two atoms per link, equivalent to i = 1, give rise to the required 
gauge symmetry. But enlarging £ here involves many atoms of the same two species 
- condensates - rather than many degrees of freedom, as in [128] . The electric energy 
term of L\ is obtained by other “gauge invariant” scattering processes, between the 
bosons of a single link. 

This simulation proposal has been treated numerically, using matrix product 
states El!- 

















CONTENTS 


37 



Figure 9. "Exact gauge invariance" - the Schwinger model example, (a) The 
desired interactions § are obtained as the F preserving boson-fermion scattering 
processes, utilizing the overlap of bosonic and fermionic Wannier functions on the 
link (b). By appropriately choosing the hyperfine levels repre sentin g t he bo sons 
and the fermions (c), only the gauge invariant interactions ( |128| ) of {TT 3 } are 
possible (d). 


6.2. SU(2) - Yang Mills theory - on-link interactions 

6.2.1. Left and right Schwinger representations. One can utilize this type of 
interactions for other gauge theories as well, such as the SU{2) lattice gauge theory. 
In m, we have suggested a way for realizing an SU( 2) lattice gauge theory, in 1 + 1 
dimensions, utilizing the prepotential method [1401114K : 142i U43| for representing the 
gauge degrees of freedom (which generalizes the Schwinger representation). In this 
method, every link is decomposed into two parts - left and right, corresponding the 
left and right degrees of freedom. The Hilbert space of each edge of the link consists 
of two harmonic oscillators, generating a separate Schwinger algebra: ap 2 on the left 
side, with the left SU( 2) algebra p4| generated by 

L c* = \Y a l ( 130 ) 

ij 


and &i ? 2 on the right side, with the right SU(2) (25) algebra generated by 


b 3 


(131) 


The left and right couples of oscillators have the same total number of excitations - 


N l = a\a\ + a\a 2 = b\b\ + b\b 2 = N R 


(132) 


(which is proportional to the total angular momentum in the Schwinger representation 
(126)), and thus, relation (601 is satisfied and the correct Hilbert space, with three 
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Figure 10. Structure of the non-Abelian lattice atomic simulator in the 
prepotential method: the simulated link is decomposed into two pieces, tied 
together by an auxiliary fermion. 


degrees of freedom and left and right SU(2) algebras, is achieved. The rotation 
matrices’ elements are constructed out of creation and annihilation operators of these 
oscillators. One can decompose them to U = UrUr, with the matrices 


U L 


1 

VW+1 


a[ -a 2 

CL<2 CLi 



1 

VNr + i 


(133) 


defined on both sides of the link. 

The implementation presented in |146| involves bosons on the links and fermions 
on the vertices, as before, but now, due to the new, exact gauge invariance method, 
the hyperfine levels of the atoms must be chosen carefully, such that only the desired 
scattering processes, which correspond to gauge invariant interaction, will be allowed 
by hyperfine angular momentum conservation. The links were not established as 
"first principle interactions", but rather "glued" out of two building blocks, using 
auxiliary, ancillary fermions, in second order perturbation theory, forming effective 
link interactions: each of the building blocks is obtained using the above prescription, 
involving transition to the auxiliary fermionic site in the middle, whose energy is too 
high, and in order to satisfy this constraint, one obtains, in second order, the full SU (2) 
link (see figure 10 ). Nevertheless, gauge invariance is not effective in this proposal 
in the sense there is no Gauss’s law constraint, and this is the achievement of this 
work. The link bosons involved contain four species playing the role of the harmonic 
oscillators of the prepotential method, from which the algebras are constructed, and 
other four "auxiliary" species, serving as "reference baths" required for the number- 
conserving scattering processes. 

However, this is not a simulation of the "full" 51/(2) Kogut-Susskind theory using 
the exact prepotential Hamiltonian; Much like some of the cQED simulations (in the 
spin gauge method), in which the electric field was truncated, here it was restricted to 
the representations j = 0, This, along with the square roots in the denominators 
of prepotential representation of the group elements (rotation matrices) - that were 
not achieved in the simulation - have generated some inaccuracies compared to the 
Kogut-Susskind theory, which result in the dynamics being accurate only in a regime 
where one may consider the gauge-matter interactions as a small perturbation (to 
sixth order in it). 


6.2.2. Truncation scheme for the Hilbert space of non-Abelian groups. One can 
overcome the problem of square roots of operators in the denominator, as well as 
the need for auxiliary bosonic species and decomposition of the links, if the truncation 
scheme of 1 (13 is utilized. In this method, we still work with a truncated theory, but 
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the truncation is done in a gauge invariant way, which does not involve decomposition 
of the link into two (left and right) pieces. 

The bosonic degrees of freedom are represented by single bosonic atoms residing 
on each link of the lattice. The state of each such atom may be one of AT = 
2 

E (2j + 1) internal levels, labeled by \jmrn ). These are generated from the 
3=0 

local atomic vacuum state |0) by the creation operators a^ m ,, with j £ 

—j < 7?r, m' < j: 

\jmm!) = a^ m , |0) (134) 

Out of these creation operators, and their conjugate annihilation operators, one 
can construct the SU (2) algebras 

( T ^)nm tt hm' R « = ( T a) m >n> a ™n> (135) 

3 3 


where a summation is assumed on double indices. Xare the j th representation 
matrices of 577(2) (for example, T a = <j a /2). These operators satisfy the left 
and right SU( 2) algebras presented above. Ladder operators may be defined as 

L ± = L x TiL y = ( T ±) «nm' and R ± = R x±i R y = E«W ( T ± ) , , C ■ 

j \ / mn j \ / m'n' 

All these operators satisfy the algebras 

The commutation relations with the creation operators are 


T 

u mm' 

7? 

' «, n nim' 


Jj 


= (n ) 7 

= *'L> ( T i) n ^ 


(136) 


implying that they undergo left and right rotations within the j representation. Thus, 
if we perform a rotation V £ SU( 2) on the left and W' £ SU(2) on the right, these 
operators will transform according to 


Jj 


V j a* j ,W^ , 

mn^nn' n'm' 


(137) 


with , IL’ll the jth matrix representation of the group elements V, . 

Next, one would like to exploit these transformation properties for the 
construction of rotation matrices - i.e., SU(2) group elements. This is done, as 
explained for general groups in m » using the Clebsch-Gordan series. We denote 
the Clebsch-Gordan coefficients by ( J,K) = (JjMm\KN) , and define 

< m ' (J, K) = < MJV (J, K) cC' N ' ( J , K) a)« N ,a J MM , (138) 


The u 3 mm , (J, K) operators undergo rotations according to the jth representation, and 
thus the corresponding rotation matrix elements [ 165] may be constructed out of them: 



2J+1 , 

- u 

2K + 1 mm 


(J,K) 


(139) 


However, as the u J mm , (J, K) operators undergo rotations independently of each other, 
we can also construct objects with the desired transformation properties without 
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- o, ±1, ±2 


m F = ± 1/2 


Figure 11. The simulating lattice structure. As usual, the bosonic potential’s 
minima are on the links, and the fermionic - on the vertices. The different shapes 
of the bosonic minima are schematically presented. 


including an infinite number of levels, i.e. with J max < oo. This is the result of 
truncating the sums in the generators and the rotation matrices, limiting them to 
J < Jmax- Thus, also for a finite J mQX one can obtain rotation matrices operating on 
the j < J max sectors, and this allows the construction of a truncated but yet gauge 
invariant lattice gauge theory, as we shall describe next. The unitarity is lost, and 
now 

tr (U^U j ) = tr (CPtPt) =2j + l- fj (J max ) P Jmax (140) 

where Pj = a mm' a mm' * s a projection operator to the j sector (recall that we are 

m,m' 

dealing with the Hilbert subspace of a single atom). 

The electric Hamiltonian will simply consist of Pj projectors with the right 
prefactors, corresponding to linear terms in the bosonic number operators. The gauge- 
matter interactions are, however, more complicated. 


6.2.3. Realization of the truncated on-link interactions. Next we shall give an 
example for the realization of the above truncation, with the j = 0,1 /2 representations. 

The scheme we were studying contains bosons belonging to a F = 2 hyperfine 
manifold, with five different hyperfine levels (as required for J max = 1). They are 
prepared such that each link is populated, as required for our representation, by a 
single atom, and the bosonic optical wells are deep enough, such that no hopping 
of bosons between different bosonic sites are possible. Thus the system remains in 
the local single-atom Hilbert subspace. The bosonic creation operators are {hj^}, for 
—2 < m = tuf < 2. The trapping potential’s minima of the different bosonic species 
have different shapes (generated by different internal potentials) - there are at least 
three different such shapes, for \m\ = 0,1,2 (see figure 111. 
to an F = 3/2 manifold, with four different hyperfine levels. 


The fermions belong 
These are required 


for the mapping between F conservation in atomic collisions and gauge invariance. 
Each vertex (fermionic site) may be occupied by two fermions at most, and in an 
alternating manner: the even sites may be occupied by mp = ±3/2 fermions, while 
the odd ones by mF = ±1/2 only. We label even vertices’ fermions by if and the odd 
by x, where m F (ipi) = -3/2, m F (if 2 ) = 3/2, m F (%i) = 1/2, and m F (x' 2 ) = -1/2. 
The fermionic sites are designed such that the fermionic wavefunctions of the two 
edges of a link will overlap on the middle on the link, also with the bosonic wave 
functions (see figure 111. 
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mr 


(b) 

Energy 

© 

H(A + A) 

9 

Vi(n 2 - A) 

-3/2 

3/2 

- %(f2i + T2j) 

© 

-h(A-A) 

F = 3/2 

© 


Figure 12. The energy levels of the (a) F = 2 bosons and the (b) F = 3/2 
fermions , as described by fl q . 


Using lasers or magnetic fields, depending on the atomic realization, the bosons 
and fermions are subject to the Hamiltonian 


H n = b ™bm + ^ + ^) 


L m^O 


+ ( x i xi - X2X2) 


(141) 


where L is the set of links, V e of even vertices and V 0 of odd ones, and 12 1 7 ^ O 2 , 
2tti 7 ^ 12 2 (see figure [ 12 } . 

Gauge invariant interactions, according to are of the form ijr^Uxi where U 
is the group element, which is a matrix of bosonic operators, as given by ( |139 1 with 
j = 1/2 and J max = 1/2. These interactions are constructed by the on-link fermion- 
boson scattering. The chosen values of tuf allow only specific scattering channels, 
thanks to the conservation of hyperfine angular momentum. By the introduction of 
Hq and energy conservation considerations, we further reduce the number of possible 
scattering processes (assume 12^2 are large enough compared to the other energy 
scales, and perform a rotating wave approximation with respect to Hq). In first 
quantization terms, the S-wave scattering is governed by a pseudo-potential Vs whose 
matrix elements are 


(142) 


2 ,m b ,-,m f ) = 


^ 2 , m' b , m' f 

7 

E C F x ^2, m' b , ^,m' f \F, m! b + m'^J 
x (^F,m b + m f \2,m b , 


(The Kronecker delta is redundant due to the Clebsch-Gordan coefficients, but is 
explicitly written to manifest angular momentum conservation.) 

By properly choosing the Cf coefficients and the overlap integrals (by 
manipulating the shape of the optical potential and hence the Wannier functions) 
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one may obtain the required link interactions as in (1391 and the electric energy as in 
equation (30 ) (which is simply a sum of bosonic number operators, using the projectors 
Pj defined above), as well as avoid extra interaction terms which conserve hyperfine 
angular momentum, but yet are not part of the desired Hamiltonian. Unfortunately, 
however, we have not been able, so far, either to satisfy all the required symmetries 
by the simulating system, or to prove it is impossible. An experimental realization 
might be possible for other atomic configurations, with less stringent symmetries (i.e., 
a different choice of atomic levels). 

For the simplest configuration described above (which could be generalized to 
more complicated atomic systems if required), with a theoretical simulating system 
with the right Hamiltonian parameters, one could proceed as follows: For even links 
we obtain the interaction Hamiltonian 


• II.c. (143) 

ij 


and for odd ones - 

+ Hc - 
ij 

with the bosonic matrix 


±( blbo + blb- 2 -blbo + blb-!} 

Vz V -bhbo + bib 1 blb 2 + bl a bo J 


(144) 


(145) 


The next step is to map between the atomic degrees of freedom b m to the ones of 
the simulated theory - a 3 mrn ,. The mapping is as follows: 


(i) For even links, 


(ii) For odd links, 


>t _ ti/2 

°±2 — a ±l/2,±l/2 

• t _ U/2 

b± 1 ~~ a ±l/2,=Fl/2 


h t 
°0 


= a 


to 

00 


h 1 

°±2 ~ 
fell = a. 


- „U/2 
“ °Tl/2,Tl/2 
11/2 


±1 — “±l/2,=Fl/2 


bl ~ a oo 


to 


(146) 


(147) 


Using this mapping, one obtains that M = U for even links, and M' = U for the odd 
ones (with U being the truncated rotation matrix, acting only within the 5-dimensional 
Hilbert space of j = 0, 

Thus, if we relabel all the fermions as \I/, we obtain the desired interaction 
Hamiltonian, 

Hint = ^2 + H.c.) (148) 

* n 




CONTENTS 


43 


6.3. Plaquette interactions: d + 1 dimensional Abelian and non-Abelian theories 


The non-effective simulations may be extended to further dimensions. The nontrivial 
generalization is from 1+1 to 2+1 dimensions, as it involves the inclusion of the 
plaquette interactions. In the previous methods, which involved effective gauge 
invariance, these were effectively obtained out of the gauge-variant interactions of the 
original Hamiltonian, using the gauge invariance constraint. Here they are obtained 
effectively as well, but out of already gauge invariant building blocks - the on-link 
gauge-matter interactions, which are not effective here, which is a great advantage. 

In [HI], we have introduced the "loop method", which uses heavy ancillary 
fermions - constrained to populate only some special vertices that may move only 
virtually, closing loops along plaquettes, carrying the gauge degrees of freedom along 
the way (utilizing the on-link interactions (128)) and thus effectively forming the 
plaquette interactions (see figure 131. Although these interactions are obtained as 
fourth-order processes and might seem too weak, they are practically only second order 
interactions, since the perturbative series of the gauge-matter Hamiltonian contains 
only even orders (fermions must hop virtually across an even distance, in order to finish 
in a vertex in which they may stay) and thus the fourth order is the second-leading 
one. 


One can introduce yet other fermionic species (non-constrained) to serve as 
"regular", dynamical fermions, on top of the ancillary static ones. 

We have utilized the loop method for a 2+1 dimensional compact QED in 
the spin-gauge approach - plaquettes are constructed from the 1+1 dimensional 
link interactions, generated due to hyperfine angular momentum conservation. This 
included a numerical proof of principle, considering some gauge invariant but non- 
Kogut-Susskind corrections to the Hamiltonian, due to the effective calculation in 
the truncated case. It was shown that the effective Hamiltonian still gives rise to the 
expected spectrum regardless of these corrections in several regimes of the Hamiltonian 
parameter, also outside the strong coupling limit. 

The loop method has also been formulated and utilized for "ideal" simulations, 
rather than the above described "realistic condition" cQED in 2+1 dimensions. This 
includes "ideal" (not truncated) cQED, as well as SU(N ) gauge theory with unitary 
matrices (which is not realistically achieved in the current simulation method). 

Finally, we shall comment on the simulation of the Z^r lattice gauge theory, whose 
relevance to QCD confinement has been stated above. Since this theory involves finite 
Hilbert spaces on the links, one could, in principle, simulate the exact model, with 
no truncations and approximations. An explicit example for Z 3 was given. Since 
the local Hilbert space is finite, we use different bosonic levels on each link, playing 
the role of the different eigenstates of the Z ^ P operator. The on-link interactions 
are with auxiliary fermions (which are not real Zjv, and are eliminated in the loop 
method), which can thus be also hard-core bosons (the statistics is irrelevant). The 
gauge invariant operators are obtained using the conservation of hyperfine angular 
momentum in collisions, as before, but the challenge here is to realize the cyclic 
"Escher staircase" nature of the group. This is done by hybridization of atomic levels, 
as discussed in [Mi- 

Dynamical fermions may be included in these simulations as well, following the 
formulation of [ 163 . 
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Figure 13. The loop method: Heavy, ancillary fermions (of two types, 
represented by the color) are constrained energetically to occupy certain vertices. 
Virtual loops of such fermions (as in the middle plaquette), carrying the gauge 
degrees of freedom along with them, generate the plaquette terms. For further 
details refer to m- 


6-4- Gauge invariance from atomic interaction symmetries - some general remarks 

To conclude, the approach of exact gauge invariance seems to be better than the 
previous effective approach, as it is applicable to a larger class of models (with the 
loop method for the plaquette interaction, for example). 

This is since it does not have to include the Gauss’s law constraints, which are hard 
to realize (and even harder for non-Abelian gauge groups, due to the non-commuting 
and complicated generators) and thus this approach seems the preferred to proceed 
with. 

Regarding the simulated theories: the simulation of cQED in the non-effective 
approach seems closer to a physical realization than the previous, effective ones. A 
better physical realization of the non-Abelian models is required, in order to 

(i) Allow for a better 1+1 simulation, with no need to use perturbation theory at all 
(although the current realization does not rely on perturbation theory for gauge 
invariance, only to the generation of links). 

(ii) Include further nonperturbative regimes of 1 + 1 dimensional theories. 

(iii) Extend the simulations to 2+1 dimensions, in a way that plaquettes will be 
generated using the loop method even out of gauge invariant building blocks 
which contain non-unitary matrices; This should be possible using the truncation 
scheme of m, but as explained above, a concrete physical realization is still 
lacking. 

(iv) Allow for simple generalizations to SU ( N ) with N > 2 and further gauge groups. 
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Figure 14. In the gedanken experiment presented in f200l . whose realization 
was described in ESI- a meson, consisting of two static quarks/charges, whose 
length is initially xq , is brought into a superposition of two lengths, using two 
internal levels of one of them (e.g. It") and |4.)) which are affected by two different 
potentials, and thus follow two different trajectories; At time t the two mesons in 
superposition are brought to the same position again, resulting again with an xq 
sized meson. After performing a Ramsey spectroscopy, within a confining phase, 
the probabilities for being in either of the two external states depend on the area 
enclosed by the trajectories A. This experiment may be used for detecting a 
confining phase, and within it - for measuring the string tension. For a more 
detailed account of this procedure, the reader may refer to [2001 . 


7. Discussion 

The class of Abelian and non-Abelian lattice gauge theories that has been considered 
at the present work constitutes the building blocks of the standard model of elementary 
particles, that currently provides us with the best known description of high- 
energy physics phenomena. The standard model has been so far extensively tested, 
and currently agrees with the experimental results from particle accelerators and 
cosmological observations, up to energy scales of ~ 10 12 eV. On the other hand, 
the physics which is used here to describe the many-body system of ultra-cold atoms 
is well understood and has been tested with high accuracy, down to energy scales of 
~ lCU 7 eV: ~ 19 orders of magnitude below the energy scales of the simulated physics! 

Furthermore, the difference between the simulated system and analog-simulator 
is not only in the energy scales: there rather exists a striking difference in the formal 
structure and the physical principles that govern the physics of high-energy phenomena 
and ultracold atoms. The atomic system at hand is non-relativistic, and seems to lack 
the sort of constrained symmetries and dynamical structure that are so essential for 
relativistic HEP models. Unlike high energy physics, atomic dynamics are number 
conserving, and manifest neither local gauge invariance, nor Lorentz invariance, which 
lead, among other things, to the familiar causal space-time structure, to charge 
conservation laws, to long-range interactions in QED and to quark confinement in 
QCD. 

So the problem at hand seems conceptually not obvious: can high-energy field 
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theories be simulated by the physics of ultra low-energy systems? 

It turns out, however and somewhat surprisingly, that there are several ways 
to manifest both Lorentz (as an effective long wavelength property) and local-gauge 
invariance. The latter can be obtained in two ways: It may arise either as a 
property of a low energy sector of the atomic system - hence as an effective symmetry, 
or, alternatively, by re-arranging the interactions such that the given low-energy 
symmetries of the atomic collisions are converted into a form which is equivalent 
to local (Abelian or non-Abelian) gauge symmetries. The first case of "effectively 
emerging" gauge symmetry seems to provide yet another framework to study models 
in which local gauge invariance is not a fundamental property of the complete theory, 
but rather emerges as a low energy phenomenon as suggested in |201] . The alternative 
"exact" mapping might, on the other hand, be more suitable for mimicking theories 
wherein local gauge invariance is a fundamental property. 

In this report, we have reviewed recent progress on quantum simulations of lattice 
gauge theories using ultracold atoms, focusing on some of the works and approaches. 
These proposals serve as proofs of principle of the possibility to simulate high energy 
physics, and even gauge theories, using ultracold atoms. 

Our current results already suggest that simple enough models, such as compact- 
QED (C7(1)), and possibly also more involved, non-Abelian SU(2) models, manifesting 
exotic QCD effects, (such as quark confinement) could, indeed, be experimentally 
studied using current table-top cold atoms experimental methods. 


7.1. Advantages of quantum simulation over other approaches 

(i) As fermions "come for free" in these atomic scenarios, one may avoid the 
"sign problem" of Grassman integration and consider regimes of finite chemical 
potential, such as the exotic color-superconductivity and quark-gluon plasma 
phases of Quantum Chromodynamics |!) ll 1001 llOll . 

(ii) Real time dynamics is possible here, rather than in the statistical Monte-Carlo 
simulations (and this applies already for the current proposals). This allows, or 
shall allow, 


(a) Observing real time dynamical phenomena, including the ones described in 
the papers: deforming and breaking of flux tubes and loops, pair creation 
etc. 

(b) Changing adiabatically the Hamiltonian parameters - finding ground states 
of systems, by switching from “trivial” regimes where the ground state is 
known (e.g. no-interactions, strong limit) to other regimes; Also probing 
for phase transitions and creating phase diagrams of theories, depending on 
the theory’s constants. This shall be highly relevant for QCD once feasible 
simulation systems are proposed and realized. 


(iii) Possibility to consider finite temperature models as well^lS, SS:, I TOO I , including 
with real time dynamics . 

(iv) Possibility to realize otherwise only gedanken experiments. For example, a 
realization of the gedanken experiment proposed in |20U] . which brings the idea 
of Ramsey spectroscopy [202] into high energy physics, suggesting a way to probe 
for the confining phase and measure the string tension related with confinement 
(see figure 141. 
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7.2. Simulated models 

As described in this report, a large class of Abelian and non-Abelian physical models 
can be simulated, including, explicitly, the lattice gauge theories of 

(i) Compact QED - Kogut Susskind and spin-gauge (truncated) models, in 1+1 
and 2+1 dimensions, with both effective and exact gauge invariance, with or 
without dynamical matter. In 1+1 dimensions, the simulation is possible with no 
perturbation theory at all. 

(ii) Ztv - so far, Z 3 has been extensively studied - 2+1 pure gauge, or with dynamical 
matter (following [165] ) - an exact quantum simulation of a non-truncated Hilbert 
space. 

(iii) SU(N) - so far SU( 2) in 1+1 dimensions, with limited applicability to more 
dimensions has been studied in detail, while other SU(N) groups still await 
extended study. 

The simulatable theories, as presented in this report, are summarized in the table 
below, where KS - Kogut-Susskind, trunc. - truncated, st.c . - Strong coupling, YM - 
Yang Mills. 



1+1 with matter 

2+1 Pure 

2+1 with matter 

1/(1) 

Full KS + trunc. 

Full KS + trunc. 

Full KS + trunc. 

z 3 

Full 

Full 

Full 

SU( 2) 

YM + trunc. 

YM + trunc. (st. c.) 

YM + trunc. (st. c.) 


For other Zjv and SU(N) groups, the simulation formalism still applies - however, 
experimental realizations seem highly challenging, and thus, perhaps, one should seek 
for other methods. 

7.3. Open problems 

These models and simulation proposals still require some more study, both from 
theoretical and experimental aspects. 

On the theoretical level, it is very important to understand the implications 
of using truncated systems - are systems with a finite local Hilbert space accurate 
enough for simulations of Kogut-Susskind models, if one exploits an appropriate 
coarse-graining scheme? And, resulting from that, what should be the continuum 
limit of such theories, if it exists? 

Such truncated models may be treated using tensor network techniques, and 
indeed, some studies approaching lattice gauge theories with tensor networks have 
already been carried out [T7210773112H!T112TR1 fT7^1 Etna KOBl 1 17B 11 17f> I H77112071IT7B1 !2nB| . 

As far as the simulating systems are considered, the following should be studied, 
in order to progress on the experimental side: 

(i) Techniques for generating complex lattice and superlattice structure - both for 
more challenging simulation models and for increasing the number of atoms which 
may populate a single site. 

(ii) Improving the controllability of parameters - new ways to control scattering 
coefficients, for example; Also, Strengthening the possible interactions is 
important, for example for the realization of effective, ~ J 2 /f7, interactions, or 
the second SU(2) realization scheme presented above, as well as, potentially, the 
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use of the truncation scheme of [163, for the simulation of other non-Abelian 
theories (with either compact Lie or finite gauge groups), including, of course 
SU{ 3) for QCD. 

(iii) Methods to create longer-living and more robust lattices - cooling techniques, 
ways to deal with atomic losses and more. 

(iv) Measurement techniques - the current single-addressing techniques [ 1891 : 19011191 
require very sophisticated experimental methods, limited to a small number of 
experimental groups. Development of such, and other techniques, shall help 
boosting the quantum simulations of complicated physical models. 

7-4- Summary 

Atomic physics and high energy physics have, indeed, what to offer to one another, 
suggesting a very interesting and exciting physical research focusing on their interface. 
Reconstructing high energy physics from atomic components offers a great lesson in 
understanding the interactions unfolded within the standard model of particle physics. 
And on the broader sense, the continuation of such works on quantum simulations of 
high energy physics, and the growth of the community working in this direction, may 
lead to useful new ways, methods and results in understanding the most fundamental 
degrees of freedom in our physical universe. 
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